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BOUSSINESQ’S PROBLEM FOR A RIGID CONE 


By A. E. H. LOVE (Oxford) 
[Received 1 March 1939] 


1. Ly ‘Boussinesq’s Problem’ an elastic solid body, which is regarded 
as being bounded by a plane and otherwise unlimited, is deformed 
by normal pressure applied to part of the plane boundary, the rest 
of the plane being free. In one form of the problem the pressed area 
and the distribution of pressure over it are given, and it is desired 
to determine the normal displacement of a point of the initially plane 
boundary; in the other form the normal displacement of a point 
within the pressed area is given, and the distribution of pressure 
over that area is to be determined. In this second form of the 
problem the elastic solid is considered to be deformed by there being 
pressed against it a perfectly rigid body of given shape. The case in 
which the surface of this rigid body is a surface of revolution, with 
its axis at right angles to the initially plane boundary and finite 
curvature at its vertex, and the case in which it is a right circular 
cylinder, with its axis at right angles to the initially plane boundary 
and a plane end at right angles to that axis, were both solved by 
Boussinesq (1). The former case includes that where the rigid body is 
spherical, the latter may be referred to as that of the ‘cylindrical die’. 

These solutions have been found to be of practical importance (2), 
and it is desired to obtain also the solution of the problem in the 
case where the surface of the rigid body is a right circular cone with 
its axis at right angles to the plane, and its vertex penetrating into 
the region initially occupied by the solid. 

2. The initially plane boundary of the elastic solid will be taken 
to be the plane z = 0, and the positive sense of the axis of z will be 
taken to be directed towards the interior of the solid. The axis of 
the cone will be taken to lie along the axis of z. Use will be made 
of a system of cylindrical coordinates p, ¢, z, the distance of a point 
from the axis of z being denoted by p. The displacement of any 
point of the solid in the direction of the axis of z will be denoted 
by w. The value of w for the point of the solid that is at the vertex 
of the cone will be deaoted by wy. If 2a is the vertical angle of the 
cone, the equation to its surface is 

2 = Wy—pcota. 
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The strained surface of the solid will fit the cone over the part 
between the vertex and a certain circular section. The radius of this 
circle will be denoted by c. Then the pressed area is given by 





2 = 0, pe, 
and the value of w on the pressed area is given by 
wW=Wy—pcota (OS p<e). 

The pressure applied to the pressed area will be denoted by p. 

3. From the general theory (3) it is known that in axially sym- 
metrical cases the displacement and stress at any point can be 
expressed in terms of a function V which satisfies the conditions: 

(i) V is harmonic everywhere except on the pressed area; 

(ii) on the pressed area V = Aw, where A is a certain constant; 

(iii) on the plane z = 0 outside the pressed area dV /éz = 0. | 

Condition (i) implies that V?V = 0 everywhere except on the 
pressed area, that V and its derivatives of the first and second orders 
with respect to p and z are continuous in the same region, and that 
V +0 at infinite distances in the order of the reciprocal of the 
distance or a higher order. In condition (ii) the constant A is given by 

Wen So, 
1—o? 
where # and o are the Young’s modulus and Poisson’s ratio of the 
material of the solid. These constants H and o will be used in this 
paper instead of Lamé’s constants A and yw which were used in the 
paper (3). We recall the relations 
E = 2u(1+0), A(1—20) = 2yo. 


The pressure p is connected with V by the equation 


eee 
p = —— lim |—}, 
27 oz 
27 2++0\ OZ 
the notation indicating that the limit is approached by allowing z to 
decrease to zero. It is seen that V is the Newtonian potential due 
to fictitious matter distributed over the pressed area with surface- 
density equal to p. As V is determined by the conditions (i), (ii), (iii), 
p also is determined by these conditions. 
4. The most obvious way of finding V’, when the boundary of the 
pressed area is a circle, and the value of w on it is a function of p, 
would seem to be to seek a series of zonal spheroidal harmonics which 
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shall be equal to this function when z = 0 and p < c, and have a 
zero z-derivative when z = 0 and p > c. The zonal spheroidal har- 
monics will be functions of elliptic coordinates » and 6, which are 
most conveniently introduced by the equation 
z+ip = csinh(n+76). 

The variables 7, 0, ¢ are a set of orthogonal curvilinear coordinates, 
and the ranges of variation of them are similar to those of spherical 
polar coordinates, viz. n is 0 on (z = 0, p < c), elsewhere » > 0, and 
7 > © when p or z > 00; @ is 0 on (p = 0, z > 0), and z on (p = 0, 
z <0), and elsewhere 0 < 6 < 7; in particular 06 = 37 on z = 0, 
p > c; ¢ can be considered as lying between —z and z. It is strictly 
correct to regard a locus ¢ = constant as an infinite half-plane 
bounded by the axis of z. 

From the theory given by Hobson (4) it appears that V must be 


of the form a 
= A Q,,( yf 
V “2 "Q(i0) tH) 


where A = sinh n, p = cos8, 
P.,, is the symbol for Legendre’s functions of the first kind, Q,, is the 
symbol for those of the second kind, and A,, is a constant such that 
> A,, P,(u) is the Legendre’s series for Aw on z = 0, p < c. To avoid 
ambiguity the functions Q,, will be understood to be defined as by 
Hobson (5). 

5. In the cone problem, w on the pressed area is given by 

w = wy—c cot asin 8, 

and the Legendre’s series for sin@ is known (6) to be given by 
sin@ =} i-|5— it S (4m41)-m—"/1 $C a) 


422 ° 2m+2| 2.4...2m 
m=2 


Hence we get 


J 
A, = A(wo—}ae cot a), A, = A= 


» $id ccota 
242 
2m—1/1.3... (2m—3))? 


* m+2| 2.4...2m f 


an — A $ar(4m = 1) ccota (m > 2 


Aon = 9 (m D1) 


In the region z = 0, p > c we have 


= = ((-)l 
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where h = |d(yn+1760)/d(z+ip)| and @= 32, 
and this is seen to be the same as 
: {1/(csinh )}(@V /eu),.~0, 
which vanishes because dP,,,,(u)/du vanishes with yp for all values of m. 
The function V has now been determined. 
6. To find the corresponding formula for p we observe that, in 


the region z = 0, p< ¢, 


lim (=) — (! | ; 
2>+0 \ Cz on 7=0 


in which 0 < 6 < 3x. The right-hand member is the same as 
{1/(e cos 8)} (EV /€A),—o. 
Hence we get 
] . .Q; (10 : 
p=—-s A, i ent? )P (u) = 4 = 4), 
a=77CpL in, Q,, (20) 
% = 
in which Q/,(¢) means {dQ,,(¢)/d¢}, and is to be evaluated for £ = 7A 
and A = 0. 

According to Hobson (7) the Q,, functions, as defined by him, 
satisfy the same recurrence equations as the P, functions. From his 
definitions it is easy to obtain the particular formulae 

Q,(tA) = —icot—A, Q, (iA) = Acot—A—1, 
and from these it is easy, by aid of the recurrence equations, to 


prove that 
Z P . 1.3...(2m—1),. 
aoe — —_ m+1 1 
Q,(20) = hom, Q.,,(00) = (—1) oa a din, 


Qi(i0) , Qi, (i0) (—1) 2.4...2m 
= 1, 40) = (—1)” : 
” 1.3... (2m—1) 


The formula for p now becomes 


A [2w,,_- fij< 4m+1 ] 
) = — +e cota)|—=+ P,,, (he) |; 
‘ | 13 Pa ja 


2rcu| 7 2m—1)(2m-+2 
or, as it may be written, 


A [/2w | 
a. | (FA cot a) +e cot a| B+ (P+ BP) + (BRT AR) + ot 


27rcu | 


1 l 
ah P,..- Pht 5.3 I 
| fra amT 9m+2 2» }] 
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In the last expression the terms in the bracket multiplying c cot « 
can be regrouped, without altering their order, as 
] a 1 | | 
; 9 Fot 28)+ 3g g eat 4h) +. 
1 

—_= 2m—1)P,,,-0+ 2mP, ae 

(@m—1)2m Fama 2m Fan} + 
which, by the recurrence equation 

nP,+(n—1)P,-» — (2n—1)uP,,_, 


can be expressed as 


l l l | 
3P,- 7P,4+-...+ 4m—1)P,,,-1+..-}. 
M\T.2 iT34 ? (2m—1)(2m) ‘ , Mama | 
Hence we get 
_ A {(2w, .- = 4m—1 P | 
p= seal | - —c cota Tere > (2m—1)2m 2m—1(H) j. 
m=1 


This expression shows that, unless 
Wy = 4c cota, 


p will become infinite at the circular boundary of the pressed area. 
As the pressure will certainly not become infinite at this circle, it 
will henceforth be assumed that c and w, are connected by this 
equation. 

The expression for p also shows that p — oo when p = 0, for the 
series diverges when » = 1. This result indicates great concentration 
of stress near the vertex of the cone, which is to be expected. It 
may be remembered that in Boussinesq’s solution for the case of 
the cylindrical die there is infinite pressure at the edge of the die. 
The conclusion to be drawn would seem to be that the pressure 
against a solid of a very rigid body, with a projecting very sharp 
edge or point, is likely to produce some permanent set in the part 
of the solid that is near to the edge or point. 


7. The formula for p is now 


A = 4m—1 
p=> cot a > ( Pom—1()- 
nN 


2a 2m—1)2m 
v=1 


The series on the right is the Legendre’s series for } log{(1+-4)/(1—)}, 
for, in the first place, this function satisfies the conditions which are 
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sufficient to secure that it can be represented by such a series (8); 
and, in the second place, if the series is 


> a, P,(h), 


n is clearly uneven, and the coefficient a, for uneven n, is given by 
1 
4 (ie 1+-p 
a, = | —log P of dp, 
onl | 3*i_., (He) Op 
“1 
where the value of the integral can be proved to be 2/{n(n+-1)}. 
We have, in fact, 
1 
] 
5 log 


© 


] 
] 


ite d { 1—p* dP,) 
1—p du \n(n+1) dy | 


du 


at 


TH p (2) du sali cae 
—p 


l—p? 1/ 1 1 \dP, 2 
| =-{——_—_+4 du = (nm uneven). 
n(n+1)2\1+p° 1—p) du n(n-+1) 


1 
' A 1 
Hence we find p = —cotalog eS. 
47 l—p 

8. The result that p can be expressed in a closed form suggests 
the possibility that V also might be expressed in a form free from 
infinite series. The natural way of finding V, as the potential due 
to surface-density p, would seem to be to evaluate by direct 
integration the potential, say V,, at a point on the axis of symmetry, 
and then seek to generalize the expression so found (9). 


If we write V,, = (A/47)cot «aK, 


, , 1+ ] 
we shall have Y= | log 2 a ~,. 2xp dp, 


1—p y(p?+2 
0 
in which p = ./(1—p?/c?) and 0 <p <1. 
It will be sufficient to consider the case where z is positive. Then 
it is convenient to write 


l-+p (1+ 


log = log 


and we have 


Ny : ¥(1 + d £.l(p2-+-22)—z = 
z= | log(1 rH) we +2")—z} dp j og i 


0 0 
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Both the integrals in this can be integrated by parts, and in each 
case the terms at the limits vanish, so we get 


7 l+p dp p 
0 


As regards the first of these we have 
e 1 


I p21 22)» (22 2__¢2,,2)_-+ 
j wet) ap = [etre du. 
J ltp dp l+p 


0 
This can be evaluated by ordinary methods, and its value is found 
to be 
2] ?_} 2] I | 2 | -1 |» I »2 2 
—zlog 2+z 08 +343 ‘pias +z—,/(z*+c?). 
\(z*-+-c*) 
We also have, again by ordinary methods, 


(z2+-c?)+2z , 


lt a = —zlog’ += +,/(22+¢?)+2 log 2—z. 


; P 

0 
Hence we find, when z > 0, 

2 B 5 +etan-1°\, 
| (2?+c?) 2) 

9. The key to the method of generalizing this formula so as to 
obtain an expression for V which shall be valid at all points in the 
region z > 0, whether they are on the axis of z or not, is found in 
the fact. that, if r denotes distance from the origin, the function 


V,, = A cot alzlog 
‘ 


zlog(z+r)—r, 
sometimes called the ‘second logarithmic potential’, is a solution of 
Laplace’s equation V? = 0. At points on the half-line given by p = 0 
and z > 0 this function becomes 
2(logz+log 2—1), 
so we may generalize the term zlogz which occurs in V, as 
zlog(z+r)—r+2(1—log 2). 
It occurred to me to try adding together three second logarithmic 


potentials, one of them the above function related to the origin, and 
the other two related in the same way to the points 


p= 9, z= +1. 
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Let R,, R, denote the ‘complex distances’ of a point (x,y,z) from 
these points, the phases (or arguments) of R, and R, being adjusted 
so that, on the half-line (p = 0, z > 0), R, and R, may be equal 
respectively to z—ic and z+ic. If we write, as a definition, 


W = {zlog(z+ kes te }{(z—ic)log(z—ic+ R,)— R,}— 


— H(z+ic)log(z+ic+ R,)— R,}, 
the value of W on the half-line will be 


z[log 2z—} log{2(z—ic)2(z+-ic)}]+ dic log <— 


or zlog —+ctan-! 


V (+ Fa 
So that we have, on the half-line, 
V = AW cota. 

10. We are going to show that the function AW cot « satisfies the 
conditions laid down for V in §3. 

We begin by explaining how the moduli of R,, R, are determined 
and how the phases of R, and R, are chosen. So far we have the 
definitions 

DS. 8.1 fe__de)8 ee 
Ri = p?+(z—10)?, R: <i en 
from which we see that, if we take R3 = |R,|*e'", we can have 
R? = |R,\%e-*. With this notation 
|R,|?cosv = p?+2?—c? = |R,|? cosv 
and R,|?sinv = 2zc = |R,|? sin v. 
Also |R,|? - 
of the point (x,y,z) Pati ay salle S, H sidash are on the focal 
circle, z = 0, p = c, S being in the axial half-plane through (2, y, z), 
and H being in the other half of the same axial plane. For we have 
2 ef? — | 21 of eo oe 
| Ry |* = (p?+-2?—c?)?+ 42%c U(p—e)?+2*}{(p+c)? +27} = rir. 

In the region z > 0, with which alone we are concerned, sinv is 
positive and cos v is positive or negative according as r > c or r < ¢, 
so that there is always a value of v between 0 and z. We choose 
this value to express the phase of R3. It is determined without 
ambiguity by the condition of lying between 0 and z, together with 


the equation 
tanv = 
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This equation can be written 
p?+(z—c cot v)? = c* cosec*v, 
and it shows that the point (x,y,z) lies on a circle, whose centre is 


at p = 0, z = ccotv, and whose radius is ccosecv. The coordinate 


z being positive, the angle v is the angle subtended at (x,y,z) by 
the line SH. 


We choose the phases of R, and R, to be —4v and $v. Then the 
values of R, and R, on the half-line p = 0, z > 0 will be z—ic and 
z+ic. In general we shall have 

R, = J(r,72)(cos v—isin 4), R, = \(r,72)(cos 4v+‘sin $v), 
so that 

R,+ R, = 2y(r,72)cos $2, R,—R, = 2iJ(r,72)sin Ww. 
These can be expressed in terms of the elliptic coordinates y, 6. For 

{2,/(r, r2)eos $v}? 

= 2r,r,(1+-cosv) = 2r,7.+ 
and 
{2,/(r, r2)sin $v}* 
= 2r,r.(1—cosv) = 2r,r.—(r?+1r}—4c*?) = 4c?—(r, 
Also r,+1r, = 2c cosh n, T.—1, = 2csin#. 
Hence we have 
R,+ R, = 2csinh n, R,— R, = 2iccos 8, 
the signs of R,+ R, and —i(R,— R,) being positive. 
The product R, R, is c2(cosh*n—sin?@) or c?(sinh*-+-cos*6). 
11. We can now express W in a real form. We have 


W = zlog(z+r)—4<log{(z—ic+ R,)(z+ice+ R,)}+ 


r+3(R,+R,). 
Also 
(z—ic+ R,)(z+ic+ R,) = 22+2(R,+R,)+ R, R,+ice(R,—R,)+c? 
= c*%(sinh?» cos?0+ 2 sinh?y cos @+-sinh*y + cos*6-+ 2 cos 0+ 1) 
= c® cosh?y(1-+-cos @)*. 


Further, 


z—ic+R, _ ,_ (sinh y—1)(1+cos @) 
log ict R, a 108 (inh »-Fi)(1-beos 0) 


—2icot—sinh , 
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the principal value of the logarithm, and the value of the inverse 
cotangent that is between 0 and }7, being chosen. Hence we get 


: 2--+r 
W = zlog as 
~ ¢ cosh 7(1-+cos 8) 


+c cot“ sinh »—r-+c sinh ». 


On the half-line p = 0, z > 0 this reduces to the value found 
before for (1/A cota«)V, for on this half-line r = z = csinhy and 
cos@ = 1. 

On the plane z = 0 in the region p < c, it reduces to 

s7C—p, 
for, in this region, z = 0, 7 = 0, r = p. 

On the plane z = 0 in the region p > c, it reduces to 

csin—(c/p)+./(p?—c?)—p 

To approximate to W when r > 0 and z > 0, we note that 


r = C,/(cosh*7 —cos*#), 


1 0820 \4 
z+r = ccosh “ ]— : ——)* cos 0+ ieince 7 ay 
| cosh? cosh?) J 


which is approximately equal to c cosh y(1-+-cos 6)(1—4 cos 6 sech?n). 


so that 


z+r 


Hence z log is approximately 


c cosh n(1+ cos @) 
—}csinh y cos*@sech*y, i.e. —4(c?/r)cos?é. 

Also ¢ cot~sinh 7 is approximately c/sinh y, i.e. c?/r. And —r+esinh 
is approximately 


—c cosh n| bans asc —{1l— 
| 2 cosh? . 


or it is approximately —4(c?/r)sin*@. It follows that 
lim (rW)—4c? = 0. 
ro 
12. Real forms for the derivatives of W with respect to z and p 
are obtained by differentiating the expression for W as a sum of 
second logarithmic potentials. We begin with oW /éz. 
We have 
ow . — : 
= = log(z+r)—}log{(z—ic+ R,)(z+ic+ R,)} 


z+r 
log : 
ec cosh n(1-+cos 6) 
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On the plane z = 0 in the region p < c, this reduces to 
log{sin @/(1-+-cos @)}, 


which is the same as —}log{(1+-)/(1—y)}. 

On the plane z = 0 in the region p > ¢, it is zero. 

The function AW cota has now been shown to satisfy the con- 
ditions laid down for V, with the exception of the condition that the 
function and its derivatives should be continuous. This, however, 
is evident from the logarithmic form of W, for each of the three 
second logarithmic potentials satisfies these continuity conditions in 
the region with which we are concerned. We can infer that 


V = AW cota 


is the solution of the cone problem. 

We have noted above that the value of @W/éz on the pressed area 
is that required for —éV /éz to be equal to 2zp. 

The resultant pressure P on the pressed area is obtained by 
integrating p over the area. We have 


oa 
p= A cot « log l TH 2mp dp 
4n 1—p 


1 
= 4Ac* cota | log aarts du. 
0 


1 


1 
| . I+p e l+p d 1—p? 
Also | og du = — | log l—p du 2 du, 
0 0 


which is found easily to be equal to 1. Hence 
P = }Ac? cota = }nc? cot af H/(1—o*)}. 

This result gives c when P is known. 

The value of W outside the pressed area shows that the profile of 
the strained bounding surface of the elastic solid is given by 

w = {csin-1(c/p)+./(p?—c?)—p}cot a. 

This result, with the relations between wy, P, and c, and the 

formulae for p and V, constitute the solution of Boussinesq’s problem 


for a rigid cone within the meaning usually attached to such solu- 
tions. 
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13. For a complete solution it would be necessary to calculate the 
components of displacement and stress at any point. Formulae for 
this purpose are known (3). To apply them we must obtain real 
expressions for some additional derivatives of W. 

From the logarithmic formula we have 


OW [2z _,\p lf z2-t _1\2 
dp \etr dr 2 R, R, 


Here 


Hence 


r—csinh 7 
; : 
The limit of this as p > 0 is zero, and that of p-\(@W/ép) is finite 
except at z 0. 
Of the second derivatives of W only two need be found, viz.: 
W /éz* and é*W /épéz, for 6?W /ép? is given by 
OW 10W , ew _ 0 
op” 7) dp a2? 
From the logarithmic formula for W we found 
ow 


= log(z+r)— 4 log(z—ic+ R,)—} log(z+ic+ R,). 
Oz t 
l sinh 7 


ar c(cosh?n —sin?@) 


] ] ] ] 
ETE Este hz) 
R,—z+ice 1 R,—z—ic 
. “a. £ 


(+z 1)_ 1/1 1] 
R,'R,) rj 2k, &, 


sinh 7 iF cos 6 
2 


c(cosh?7 —sin?6) cosh? — cos?6 


cosh*ycos@ 2 
cosh?n—sin?@ r}" 
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14. The component displacements U and w in the directions of 
increase of p and z are given by 
U = 159/595)" 4-5 
2rE | ép 
l+of, >. oF) 
— {2(1—oa)V—z— }, 
27 E | ( ) éz| 
where @x/ép is the function determined by the conditions 
o|[ ex oV 
a) —|p—) = —p—, 
ss alee) < 


(6) a(e2) =p 


oz\ Op 


eV) 
C, 


ep)’ 


16 2 


(c) in the region (z = 0,p >), px = P. 


From these conditions it is found that 


6 
re | (-°S ep rant a de, 


A cota dz 00? tp 2 


ia 
» being constant in the integration. After some reduction the result 
is obtained in the form 


(: —sinh »)sinh 7 cos 6-+-cos 8 + 
c 


+ cosh? sin*@ log — en =i) 
The stress-components that do not vanish are given by 

a eV) 

ial. * ep?) 


be 
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The only derivative in these for which no expression has been given 
is 6*y/ép?, and this can be found from the condition (a) above. 

We have thus obtained formulae from which all the components 
of displacement and stress at any point can be calculated if desired. 
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SUMMARY 


When a rigid cone is pressed point foremost against the surface, 
plane when unstrained, of an elastic solid, the solid is held strained, 
as if by pressure applied over a certain circular patch, the ‘pressed 
area’. The normal displacement w of a point on this area is deter- 
mined by its value wy, at the centre and the condition of fitting the 
cone. The pressure p at a point of the patch is to be found from 
the condition that the potential V due to surface-density p on the 
patch is a known multiple of w on the patch. From this condition 
V is determined in terms of w, and c, the radius of the patch, as an 
infinite series of spheroidal harmonics. The corresponding formula 
for p is deduced, and the condition that p must not be infinite at 
the boundary of the patch leads to a relation connecting w, and c. 
The series expressing p is summed, a simple closed form being 
obtained. From this form for p a formula for V at any point on the 
prolongation of the axis of the cone into the solid is deduced. This 
also is a simple closed form. It is generalized so as to give a closed 
form for V anywhere in the solid, this form being the sum of three 
‘second logarithmic potentials’, two of them containing imaginaries. 
This is then expressed as a real closed form, and it is verified that 
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this form actually does represent the potential due to surface-density 
p on the patch. By means of the closed form of p the resultant 
pressure P between the cone and the solid is expressed in terms 
of c, so that c is known when P is known. By means of the closed 
form of V the shape of the part of the deformed surface that is out- 
side the pressed area is obtained. To complete the solution, formulae 
in terms of derivatives of V, all expressed in real closed forms, 
are given for the components of displacement and stress at any 
point in the solid, so that these components can be computed if 
desired. 








TRANSFORMATION OF HYPERGEOMETRIC 
INTEGRALS BY MEANS OF FRACTIONAL 
INTEGRATION BY PARTS 


By A. ERDELYI (Edinburgh) 
[Received 9 March 1939] 
1. THE transformation into each other of Euler’s hypergeometric 
integrals representing F eae ), viz. 


I'(y) 
B)l'(y—B) 


f. pB-1(1 —x)”-B-1(1—arz)-* dx (1) 


and spe | x*-1(] —x)¥-*-1(] —az)-F dx, (2) 
D(a) P(y—a) 
0 
has been the subject of investigations by several authors.* In a 
recent note+ I proved that this transformation can be effected by 
means of fractional integration by parts. Moreover, fractional 
integration by parts can also be used in order to obtain a more 
general transformation of (1), yielding the functional equation 
ry) f 
7 er a A-1(] —_»\y-A-1 ; 
F(a, B;y;z) = Pwr) . x\-l(]—2) F(a, B;A;az) dx (3) 
’ . 
[R(y) > RA) > 0; z Al; jarg(1—z)| < a] 
of the hypergeometric function. 
In this article I propose (i) to derive (1), (2), (3) directly from the 


definition ' 
(x),(B) 
(a, Biy52) = > ree (4) 
2 a 
of the hypergeometric function} using the conception of fractional 


* B. Riemann, Gesammelte mathematische Werke und wissenschaftlicher 
Nachlass (Leipzig, 1876), 62-78 ; Schellenberg, Dissertation (Géttingen, 1892) ; 
W. Wirtinger, Sitzwngsber. Akad. Wiss. Wien, 111 (1902), 894-900; E. G. C. 
Poole, Quart. J. of Math. (Oxford), 9 (1938), 230-3; A. Erdélyi, ibid. 8 (1937), 
200-13 and 267-77. 

+ See above, 129-34. 


+ (oj, = 1, (x) = ~ a+l1)..(a+n—1) (n 


The index of summation runs in all series from 0 to oo. 
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derivatives, and (ii), by transforming (1), (2), (3) by means of frac- 
tional integration by parts, to obtain some functional equations and 
integral representations of F. 

The methods used here also apply to the corresponding problems 
in the theory of generalized hypergeometric series.* For the sake of 
simplicity, however, I shall confine myself in this article in general 
to the ordinary hypergeometric function. Some limiting cases of the 
formulae obtained in this article yield integral representations for 
the confluent hypergeometric function. 


2. We write the rule for fractional integration by parts} in the 


form b 


b 
r d’y 3 du 
' he ane , lx. ~ 
| : d(b—2)’ sis | : d(e—a)’ ‘ : (5) 


a a 


The fractional derivatives occurring in this rule can be defined by 
integrals, if the real part of v is negative. Thus 


zx 
. 


d’u l 


= — —v-l 
d(x—ay I(—»), (w—y)“tu(y) dy 
= [R(v) < O}. (6) 
cs = . j —y)-¥-ly 
d(b—x)” as T(—») | (y—2) v(y) dy 


x 





If wu and v are expressible by means of series of the types 
u = > A(x—a)Pr—t, v= > B(b—x)°*—, (7) 


then the fractional derivatives are obtainable by differentiating these 
series term by term and using the definition 


ut) (jut 


dw™ — T(u—) (8) 


for fractional derivatives, which holds for all values of v excepting 
vy =p. Obviously (6) and (8) are in accordance in those cases in 
which both definitions have a meaning. 


* e.g. equations (1.5), (2.6), (3.2), (5.2) of my paper Quart. J. of Math. 
(Oxford), 8 (1937), 267-77 are obtainable by the method of § 3 of this 
article. 

+ E. R. Love and L. C. Young, Proc. London Math. Soc. (2) 44 (1938), 
1-34. 

3695.10 N 














178 A. ERDELYI 


3. Now we can derive (1) by means of fractional derivatives. 
Using the definition (8) we obtain 











> (x),(B),, _ T'(y) dB-y | i. 
Bla, Bsy32) = bs (yr! —«sT(B)” wilt "TF 
DO) 1 © pa (1—z)-*}. (9) 


Tip) der" 


Now, using for the fractional derivative on the right of the last 
equation the definition (6), we at once arrive at 


I'(y)z-” 


be ae 
PA —A) J ” 1(1—y)-*(z—y)?-B dy. 


F(a, B; 32) 


Introducing a new variable of integration 2 by the equation y = zz, 
the last relation turns into (1). A very similar analysis yields (2). 
Instead of (9) we can transform (4) as follows: 


ee ee — (x),(B), ——— = r (y) d- Y{ ),(B),. gAtr- | 
Fa Biyi2) = > rl * Ta) a “2% a 
— »zl-y l'(y) d\-y {yA- 
~ TA) daw TP (a, B5A;2)}. (10) 


Using again definition (6) of fractional derivatives we obtain 


vcs at en 
(P5752) = TOPO) . 


0 


(z—y)”*-ly-1F (a, B; A; y) dy. 


This equation yields, on putting y = z2, the relation (3). 
The formulae (9) and (10) are generalizations of well-known 
formulae for positive integer values of B—y and A—y respectively. 
4. Transforming (1) by means of fractional integration by parts, 
I shall now prove the formula 


1 
7 ee. l'(y) , A=—1/ 3 mi\¥—A—2/ 1 ___ are ed 
F(a, B34 32) —— P(A) P'(y—a) | x (1 x) (1 x2) 
0 
—p ' il > D (l—2)z 
< F(a—a’, B;A; xz) F (a ,B—A;y—A; — dx, (11) 


which is valid provided that 


R(y) > RA) >O and zl, |arg(l1—z)| <=. 
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Using the definition (8) of fractional derivatives we see that 
(1—a)v-B-1(1 —az)-« 


= (1—2)-*(1—ayp-F4( 14 =) 
1—z 


-aVy—B) BA fm (a')ly—B) 2% Vp vy rae 
T'(y—A) d(1—a)P-| (y—A),r! (0 x)” : 


— a—zy-« EB) FO fy yy va ply’ Be y—2: U2 
=“ I(y—aA) da—xp\" said Fy Bsy—d; z—1 : 





= (1-2) 








Now* F(a,b;c;w) = (1—w)-*F (a, 0—b5e; = 
w— 
and therefore 
(1—a)”-B-1(1—az)-“ 
y—p) BA | ae Ewe (1—z)\) 
= 1—x)Y\-(1—az)-*' F(x’, B—A; yp—A; "II. 
I'(y—aA) da—axP\' se ii “P lie =|} 
Using this result we can write 
1 
. I'(y) [ 
F(a, B;y32) = _ |. eB-1(1 —a)y-B-1(1 —arz)-* dx 
om PG y—B) J 
ry) f 
= Y [ 28 1—az)*-*{(1—ax)v-P-1(1—arz)-*} dx 
rere)” | ita 


1 
ie I'(y) oB—-1/ ] — r>)a’—a y 
= T@ly—a [ « (1—zz)* -* x 
i 
X 11 —2)P | 


Integrating fractionally by parts this becomes, by the rule (5), 


0 


(1—2)v-A-1(1 —az)- “F(x B—A;y—A; c =) dx. 


—zz || 


1 
1 l'(y) [ = ity 
»B3y32) = 1—z 1(1—2xz)-* X 
te ae 
0 
ee f 
x F(x, B—A;y—A; f =| ote M1 —a2)-9} dx. 


* E. T. Whittaker and G. N. Watson, Modern Analysis (Cambridge, 1920), 
§ 14.4, 
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Now 
dh | ? 8-1 a yaa} = (qB-A he (a—a')y 847-157] 
daB-A daB-A| Zu r! | 


= ¥ (a—a’), T(B+r) Dir-lee 
oa r! T(A+r) 7 
YR 
= y-1 ro F(a—«’, B; A; xz), 
and therefore 


F(a, B; 32) = 


I'(y) F - ve me ae 
x’ 1—a)v-4-1(1 —az)-* x 
FO —y J Oda ae) 


0 


< F(a—a’, B;A; v2) F( x’, B—A;y—A; =) dx. 


1—zz 
Thus far the proof given in this section only holds provided that 
Ry) > R(B) > 0, Rly) > RA) > 0; |z| <1, Riz) <F. 
By the theory of analytic continuation the result is true, however, 
with only the restrictions given at the beginning of this section. 

5. A few particular and limiting cases of (11) may be written out 
fully. Omitting those which yield formulae equivalent, by the trans- 
formation theory of hypergeometric series, to (1), (2), or (3), let us 
put in (11) Oa fm 


Then (11) becomes 


oo 





I'(y) ‘ ci a os 
F(a, B3 32 a%-1(] —ax)¥-*-1(] —arz)-F x 
(a, P; v3 2) Pa) (y— »), (1—2) (1—2z) 
F(«a—B, B; 4302) F(B—a, Bs ya; a) & (12) 
\ —xz 
[R(y) > R(a) > 0; 2 #1; \arg(1—z)| < z]. 


Examining the limiting case « — 00 of (11), we must make use of 


the known relation* 


ld z . 8), 2 
lim F(a. Bs75=] = ,F(6;y;2z) = >. a ah (13) 


* Modern Analysis, § 16.1. 








ee 
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We replace z by z/a in (11) and make « tend to 0. Obviously, we 


have a 
lim {1——] = 1 
a x Qa 


uniformly in « (0 < x <1), and therefore the limiting form of (11) 


“2 


runs 1 


[ 21-246; A522) dx (14) 


[R(y) > RA) > O}. 
This relation is a particular case of a similar formula for generalized 
hypergeometric series,* and also a particular case of equation (19) 
of this article. 

Deriving (14), we performed the limiting process a — 00 under the 
sign of integration. This is permissible by reason of the uniform 
convergence of this limiting process in the range of integration. 

6. I go on to transform (3) by means of fractional integration by 
parts. To do so, we first transform the right of (3) by means of the 
transformation formulat 

F(a,b;c;w) = (l—w)*-*" F(ce—a, c—b; c; w) 
of the hypergeometric series. Hence we get 


F(x, B; 32) 


—x)v-A-1(1 —az)\-9-B F(A—a, A—B; A;: 


isis y « \}\e mn) — (A—a),(A—8), A+r—ler 
1 P(A—a, AB; As 22) = ZZ Arad)? 


du-A so ),(A—B), F attr “| 


~ data Ly oy 


de fop-l —— 
— dee-\\TQ A—2,A— Pie 


* A. Erdélyi, Quart. J. of Math. (Oxford), 8 (1937), 267-77, equation (5.2). 
+ Modern Analysis, § 14.4. 
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and therefore we can write, instead of (15), 


l'(y) 
F(a, B;y3;z) = (1—ax)”-A-1(1 —az)\-9-B x 
(Bi7i2) = poe =k ny —a2p-0-B x 
ee i L 
* dat Ae (A—a, B; w3% 2) 6 : 
1 
= My) | x+-1 F'(A—a, A—B; 3 22) 
Mu). 
0 


b du-A {( —zx)vy-A-1 
* d(1—aye- | T(y—A) 
using the rule of fractional integration by parts. 
The fractional derivative occurring in (16) we compute according 
to the definition (8). We have 


(1—a)r-A-1 


(1—az)-2-#) dx, (16) 


: (1—az)-2-B 
I 4 —A 
Y , (i- x)y A-1 -_. A-a-B 

= (l—z)-=-F 7 a ) 

ie l—z 
wale z)r- a—B | z \ 1—2)v-A+-1 
rae To we i) pyre, 
and therefore 

de ((1—zyr"4 —az)\-«-B| 


d(1—ax)#-\|_ T'(y—A) 


— (1—z)\-2-B aaa Ay FE 6 Vg es 
(1 >= jor Or Lia he 


nahi . 1—x)z 
 (1—2)-a- (4-2) ee a ) 
I'(y—p) z—l 
But this is equal to 
_— y-p-1 i — 
(l—az) 2B | r) F(a+B—AA—piy—miG =) 
l'(y—p) 1—zz 
by the transformation formula of the hypergeometric series quoted 
in §4. Putting this expression into (16) we at length arrive at the 
formula 1 


F(a, B;y3;2) = 








| 
| 
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valid for 
Ry) > R(u) > 0; zl, |arg(1—z)| < z. 
Obviously (17) is a generalization of (3). 

7. Let us examine again a few particular and limiting forms of 
(17). There are many particular cases yielding (3) or formulae which 
after some transformations of hypergeometric series turn into (3). 
Omitting all such cases, let us put A= 8. Thus the first of the two 
hypergeometric series on the right of (17) reduces to 1, leaving only 

ro) 

F(«, B;y32) = y | xb-1(] —a)v-#-1(1 —xz)-* x 

Me) y—H) J 
1—2a)z 

X Fla, B—p;y— .{ 

x ( B—Bsy—-BK 
valid under the same restrictions as (17). This is a generalization of 
Euler’s integral (1), reducing to (1) if » =f. We get a similar 
formula from A = a. 

Now let us replace z by z/f in (17) and then make £ tend to o. 
Doing so we use (13) and remark that 


im (1) = ], 
p>»\ 8B 


) ae, (18) 


xz\\-a-B 
but lim (1-5) = e*, 
pa B 


Having regard to this, we obtain the limiting form of (17): 


- , 

I'(y) | 

P(n)P(y— ; 
(1) P(y—p) . 

x FF fA—a; p; —az} FF fA—p; y—p; (1—2)z} da. 


Fi x3 2) — pel l —x)Y- ler x 


Here we use Kummer’s transformation formula* 
e” ,F,(a;c; —w) = ,F,(c—a;c; w), 
which is the limiting form of the transformation formula quoted 
in §4, and get the functional equation 


1 
ss (y) j 
P(u)P(y—p) ¥ 
( 
< F, {a—A+p; pw; az} FF fA—ps y—p; (1—a)z} dx, (19) 


* Modern Analysis, § 16.11 (I). 


(a; 32) = xe-1(] —x)v-#-1 x 
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valid for R(y) > R(u) > O. 


The formula (19) is another form of the ‘addition formula’ for 
Whittaker’s confluent hypergeometric function.* 

It is perhaps not superfluous to notice that a great number of 
formulae are particular cases of (19). The simplest among these is 
Euler’s integral of the first kind, 

1 SP ae 

( aaa (| —zx)Y-# -1 dx sa v Mw) 

(y) 

0 
which results from (19) if we put «a = 0,A = yp. Among the others 
the best known are: Doetsch’s ‘Faltungs’-formulae} concerning 
Hermite polynomials, Tricomi’s addition formulat of Laguerre poly- 
nomials, Kapteyn integrals§$ involving Bessel functions, and some 
integrals|| involving parabolic cylinder functions and Bateman’s 
k,,-functions. 


8. Let us again start with (3) and remark that 


R 
a-1 F(x, B; A; xz) = > ()-(P > ler 


(A),.r! 


/ 


— d”-# & (cx),(B), D(u+r) pv tr—Lyr | 
dx BIZ, (A),. r! T(v+7) : | 


d’-- {1(u) 192 | 
= x1 F(a, B, w3A,v;%2)}. 
aT)” (a, 8, w3A,v;2 ) 


Thus we may write down (3) in the form 


C(y)P(u) 
DA)P(y—ayP(v) . 


0 


aA #(1—ax)Y A-1,, 


d’—# - “1 
> 5 oer fx”, F(a, B, 3A, v;0z)} da. 


* A. Erdélyi, Math. Zeits. 42 (1936), 125-43, equation (3.1). 

+ G. Doetsch, Math. Zeits. 32 (1930), 587-99. 

t F. Tricomi, Rendiconti dei Lincei (6) 21 (1935), 332-5. 

§ See, e.g., E. T. Copson, Proc. London Math. Soc. (2) 33 (1932), 145- 
53, § 4. 
e.g. N. G. Shabde, J. of Indian Math. Soc. New series 3 (1938), 146-51, § 2; 
N. A. Shastri, ibid. 155-63, § 5. 
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Integrating fractionally by parts, this becomes 


1 
w Brit) = TOEW Farts Bus dose) 
F(a, B; 32) = ParG—aro). x”-1, F(a, B, 3A, v; xz) X 


0 


d’- be 
haces Don yv-A- 1 
xX - qa—zy et -#(1—z) } dx. 


Now 


d¥ be 
Sy A —p (l1—2x) ~A-1 
d(1—x)-#™ 


dv’ —p 


~ d(1—2x)- = {(1—a)r [1 —(1—2) Pp} 


de { a (u—A), y—A+r— | 
a a n 1—2)? r—l 

d(1—a)-¥ | 2 r! ila 
2 ['(y—aA) 


= 1—a)v-A+H-¥-1 F(u—A, y—A; yp—A+p—v; 1—2). 
a er v) (u ? ‘i re—v x) 


Hence 

= — ‘ My)P(u) x” 1] —a)y »—+p—v—1. 
P(A) P (vy) P(y—A+p—rv) , 
< F(u—a, y—A; y—A+p—v; 1—2z)3Fy(a, B, m3 A,v;2z) dx, (20) 


and this relation is valid, by the i of petite: continuation, 
provided that 


RA) > 0, Rv) > 0, R(y—A+p—vr) > 0; z $1, |arg(1—z)| < =. 


Transform the ordinary hypergeometric series on the right of (20) 
by means of the transformation formula quoted in §4, and (20) be- 
comes 1 
F(x, B;y32) = —— Py)Ee) [ a y-1(1—a)yAtH-¥-1 x 

PA)P@)P(y—A+e—») . 


0 


x (yr y—v;3 y—A+p—y; 1— *) Fs x, B,m3A,v;~z) dx, (21) 
x 


(A) > 0, Rv) > 0, R(y—A+p—v) > 0; 2 $1, |arg(1—z)| < a], 
thus exhibiting the symmetry in A and v. 


9. In this section I write out some particular cases of (20) 
and (21). 
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Let us put v = a in (20). This yields 


1 
> I(y)P(p) Fae eee 
F(a, B;y3z) = 2-1] —x)v-4-AtH-1 y 
(a B;¥; ) hiatairtie—eln—)) | r ( v) Le 


0 
< F(u—aA, y—A; y—a—A+p; 1—2x)F(B, 3A; xz) dx, (22) 
[R(a) > 0, R(A) > 0, R(y—a—A+p) > 0; z ~1, larg(1—z)| < a]. 
We get a similar formula from v = P. 
Again, putting v = y in (21), we arrive at 
P() 
‘(u ee 1p 
»B39732) = x-1(] —a)# 4-1 F(a, B, wy, A; 22) dx 
(x, B; y32) rO)r(u—a }? ! eh" F(x, B, 3 y,A; xz) de, 
° (23) 
[R(u) > R(A) > 0; z £1; |arg(1—z)| < a]. 
This is a particular case of a functional equation of .F,.* 
A last example of this kind we obtain by putting A =: « and vy = B 
in (21). Thus, : 
(x, B34 . 2 ———— i | a%+B-y-1(] —x)v-2-B+u-1 x 
D(a) P(B)P(y—a—B+p) . 
0 





— x, y—B; y—a—B+p; Je dx (24) 
x 


R(x) > 0; R(B) > 0; R(y—a—B+p) > 0; z £1; |arg(1—z)| < a]. 


Giving suitable values to the arbitrary parameter y in this relation 





we obtain Euler’s hypergeometric integrals (1) and (2) as particular 
cases of (24). Thus, this relation is another expedient for connecting 
(1) and (2).7 

10. In order to obtain certain limiting forms of (21) we must 
deal with 


lim F(y—X.y—viy—A rey; 1-4). 


paar 
If at least one of the numbers y—A, y—v happens to be a negative 
integer (or zero), then, in the usual manner, this limit is readily 
found to be 


lim Fly—a y—v;y—A+p—p; 1— “) — Fly—Ay—v — ?) é (25) 


po 
both the hypergeometric series occurring in (25) being terminating 
ones. 

* See the first footnote on p. 181. 

T See also § 1 of the first of my articles quoted on p. 176. 








—— 
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lf neither of the numbers y—A, y—v is a negative integer, then 

(25) seems to have no significance at all, the series on the right-hand 

side of (25) being divergent for all finite values of ¢. I propose, how- 

ever, to show that (25) is valid for all values of y, A, v, if only the 
divergiah of, is defined by the equation 


F(a. -3) i POO) wa, F,(a;a—b-+1; w) + 
w 


I'(b) 
T'(a—b) 
+ T@ 


I desire to add to this definition some remarks: 


w? F,(b;b—a+1;w). (26) 


(i) If one of the numbers a, b happens to be a negative integer, 
then one of the two terms on the right of (26) vanishes, the other 
being easily transformable into the left of (26). If, for example, we 
assume a to be a negative integer —n (say), the second term on 
the right of (26) vanishes because 1/['(—n) = 0, and the first term 
can be transformed thus: 


1h L 
Tin w-" F,(—n; —n—b+1; w) 


abt (—n),27-” T'(b+n)C(1—b—n)n! (—/r 
7 ( 


T'(b) —n—b+1),r! ~ l(b) —b—n-r)r! (n—r)! , 








=>; (—n) . 1 
pon n— ( = = ol —n,b; —- 5 
(n—r) Ch r ~ 


~ 


Hence (26) is in cacaaeet with the familiar definition of ,/. 

(ii) When neither a nor b happens to be a negative integer or zero, 
then the left of (26) is known to be the asymptotic expansion of the 
right-hand side of this equation, for large values of |v}. 

(iii) Definition (26) holds only for non-integer values of a—b. For 
integer values the limiting form of (26) is to be taken. 

(iv) The reader familiar with the theory of confluent hypergeo- 
metric functions will notice the close relationship between our defini- 
tion of ,f, and the definition of Whittaker’s function W, ,,(z).* 

Now, after these considerations, it is easy to prove that, with the 
’ definition (26), relation (25) is true for all values of y, A, v. It is 
sufficient to prove (25) only for non-integer values of A—v. The 
validity of this equation for integer values of A—v then follows at 
once, both members of (25) being continuous functions of A. To 
prove (25) for non-integer values of A—v we transform the hyper- 

* Modern Analysis, § 16.12. 
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geometric series on the left of (25) by means of the transformation 


formula* 


Bes (—wy*F (a a—c+l;a—b+1; } — 


w 
1e\Ma—l ’ 
Pon (a ~ ) —w)oF /6—c4+- 1-6 ; 
' P(e—b)P (a) ; “Ww 
thus arriving at 


Fy— y—v;y—A+p—y; 1— t) 


ee P(y—A+p—v)T(A—v)/ t \” Fly eee ee t J+ 
P(u—v)l(y—v) p—t ; ‘ ' 1 - "t—p 


4 Py—At+e—v)P—A) (_t Fy vl =pdk 1p eee 
D(u—A)P(y—A)  \w—t eee 5 


Making p tend to oo we remark that, by Stirling’s formula,t 


= ety —A+e—9) 2 
-o T(u—v)(u—t)Y A 
and, by (13), 


lim Fly A, l—p+v; 1—A+»; ] = ,F,(y—A; 1—A+7;1?). 


pon t—p 
Similar relations hold, of course, with A and v interchanged. 


Hence the limiting form of (27) becomes 


t 


lim F (y- A, y—v3y—A+p—-; 1- 7} 


po 
r(A—v) 
r(y—v) 


/ 


tY \ F(y—A; l 


Ply y—Vv;— 
according to the definition (26). Hence (25) is established. 


11. I conclude by pointing out some limiting forms : fog ). To 
obtain a first example of this kind, let us replace z in (21) by z/f 
and make f tend to 00. Since 


lim F(a, 8, widivs = oM,(a, 3A, v; xz), (28) 
poo 


B 


* Modern Analysis, § 14.51. + Ibid., § 12.33. 











ON HYPERGEOMETRIC INTEGRALS. 189 

we at once get 1 

(«3 732) = ~ Dy) —— aAtv—-y—1(] —a)y-AtH-v-1 x 
P(A)P(v)P(y—A+p—») . 


) 


Fly—a y—v;y—A+p—v; 1— oFy(a, m3A,v;22z) dx, (29) 
x 


valid for RA) > 0; Riv) > 0; R(y—A+p—v) > 0. 
This integral was established in a different manner by Meijer.* 

We replace x in (21) by t/u and make yp tend to 0. The range 
of integration in ¢ runs from 0 to »; thus we obtain in the limit an 
infinite integral. We remark that, by Stirling’s formula, 

. P(y—A+p—v 
lim (y rE i, 
po T(p)py-' 


; t y—A+p—v-1 
and that lim (1-<) =e, 
B 


pore 
and recall (25) and (28). Now we can work out the limiting process, 
if we suppose R(z) <1 in order to obtain a convergent integral. 
With this restriction the limiting form of (21) is 


F(a, B; y; 2) 


= ly) [ rte lily —A ys -3) oF, (a, B;A,v;tz) dt, (30) 


D(A) P(r) 
0 
valid for RA) > 0; Riv) > 0; R(z) <1. 
Some special cases of (30) are worth mentioning. Let us put v = y 
in (30). Then we obtain 


F(a,B;7;2) = Fay | Plot, F(x, Bs y, A; tz) de (31) 
; 


[R(A) > 0; R(z) < 1}. 
Again, put A = a, v = f in (30). These particular values yield Gold- 
stein’s integralt 
F(a, B;y;z) = Ranh | (218-1Ne OR y—a y—B; -;) dt 
. (32) 
[R(a) >.0; R(B) > 0; R(z) < 1]. 
Other particular cases of the formulae are left to the reader. 


* C. 8. Meijer, Proc. Akad. Wet. Amsterdam, 41 (1938), 1113. 
+ S. Goldstein, Proc. London Math. Soc. (2) 34 (1932), 103-25, equation (54). 








THE ‘EASIER’ WARING PROBLEM 
By W. H. J. FUCHS (Aberdeen) and E. M. WRIGHT (Aberdeen) 
[Received 7 February 1939] 

1. In a former paper} Wright defined v(k) as the least value of 
such that every integer m can be expressed in the form 

n = kt aok...tak—ak,...—ak, (1.1) 
where 7, Xj, 2,..., % are integers and 0 <r <s.{ He proved that 
v(2) = 3; 4 < r(3) < 5; 8 < v(4) < 12; and 5 < (5) < 10, the re- 
sults for v(4) and v(5) being new. Our object here is to carry these 


investigations further and, in particular, to find upper and lower 
bounds for v(k) for 6 << k < 20. 


We take k > 3 and we write A(k, m,n) for the least value of s such 


that the congruence 


ak tak... tak—ak.....—xrk =n (modm) 


is soluble for some r = r(n), and 
A(k,m) = max A(k, m,n). 
n 
In the paper mentioned Wright showed that the existence of an 
identity in x 
j j 

> (x+a,)*¥— > (x+a;)* = Cx+D, (1.3) 

i=1 i=1 
in which aj, @,..., @;, Aj, @g,..., @ are integers and C + 0, implies the 
inequality v(k) < 25+ A(k, C). (1.4) 


On the other hand, a necessary condition for the solubility of 
(1.1) is the solubility of the congruence (1.2) for all m. Hence, if 
ee ee 
We WINE, A(k) = max A(k, m), 


so that A(k) is the least value of s such that (1.2) has a solution for 
all m and m and for suitable r = r(n,m), we have 
v(k) > A(k). (1.5) 
+ E. M. Wright, ‘An easier Waring’s problem’: J. of London Math. Soc. 
9 (1934), 267-72. 


t Here r depends on n. 
§ Wright, loc. cit., uses A(k) with a different meaning. 
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In $2 we develop in some detail the theory of the congruence 
(1.2) and determine various upper bounds for A(é) and similar auxi- 
liary functions. We apply these results to determine a formula for 
A(k) in certain cases, notably when A(k) > k, and to calculate A(k) 
for all k < 36. Much of our work follows closely that of Hardy and 
Littlewood on ['(k),+ though we encounter a complication which does 
not arise in their problem. We have adopted the notation of P.N. 8 
as far as possible and we omit proofs where these differ from the 
corresponding work in P.N. 4 or P.N. 8 by merely trivial modi- 
fications. 

In §3 we apply our results to determine upper and lower bounds 
for v(k) for values of k up to k < 20. For this purpose we require 
identities of the type (1.3) and we show how these may be found. 

For large values of k the best known upper bound for v(k) is given 
by Vinogradoff’s well-known result for G(/) in the ordinary Waring’s 
Problem combined with the trivial inequality 


v(k) < G(k)+1. 


2. The numbers A(k,m) and A(k) 


2.1. We use a|b to denote that a is a divisor of b, and a/b to 
denote the contrary; we write (a,b) for the greatest common factor 


of a and b. 

We write ['(k,m,n) for the least value of s such that (1.2) has 
a solution with r = s and A(k, m,n) for the least value of s such that 
(1.2) has a solution for every r (0 <r <3). We write also$ 


['(k,m) = max I (k, m,n), I'(k) = max T(k, m), 
A(k,m) = max A(k, m,n), A(k) = max A(k, m). 


n m 
Hardy and Littlewood (P.N. 8) gave methods to calculate ['(k, m) 
and I'(k), and we give corresponding methods to calculate A(k, m) and 


+ G. H. Hardy and J. E. Littlewood, ‘Some problems of ‘‘ Partitio Numero- 
rum”: (IV) The singular series in Waring’s problem and the value of the 
number G(k)’, Math. Zeits. 12 (1922), 161-88; ‘(VIII) The number I(k) in 
Waring’s problem’, Proc. London Math. Soc. (2), 28 (1928), 518-42. We refer 
to these papers as P.N. 4 and P.N. 8 respectively. 

t See, for example, Landau, Uber einige neuere Ergebnisse der additiven 
Zahlentheorie (Cambridge, 1937), Kap. 1. Vinogradoff has recently improved 
his result for k > 800 (Trav. Inst. Math. Tbilissi, 5 (1938), 167-80). 

§ This definition of ['(k) is equivalent to that of Hardy and Littlewood 
except when k = 4 (cf. P.N. 8, p. 525). 
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A(k). The work is greatly simplified by the fact that we need only 
consider the case when m is a power of a prime. This is a consequence 
of the two following lemmas. 

Lemma 1. Jf (m,,m.) = 1, and if, for given n, s, r, there is a solution 
of (1.2) when m = m, and a solution when m = mg, then there is 
a solution when m = m, Mg. 

Let 2; = a; and x; = b; (i = 1, 2...., s) be the solutions for m = m, 
and m = m, respectively. Since (m,,m,) = 1, we can find two num- 
bers d,, d, such that 

d,ms, = 1 (modmy), d,m, = 1 (mod mz). 
Then x; = d,m,a;+d,m,b; is a solution of (1.2) for m = m,mzg. 

From Lemma 1 we deduce that, if (m,,m,.) = 1, then 

A(k, my m,n) = max{A(k, m,,n),A(k, ms, n)} ) 


A(k,m mg) = max{A(k,m,), A(k, ms)} J 
and similar results for [(k,m,n) and I'(k,m). Hence 


A(k) = max A(k, p'), C(k) = max (hk, p'), 
pt p,t 

where p runs through all primes and ¢ through all positive integers. 
Nothing similar is true in general for A(k,m,n) and A(k,m); for, in 
the determination of A(k,m,n), 7 is a function of m as well as of n, 
so that Lemma 1 does not lead to results corresponding to (2.11). 
We can only say that 

A(k, m, m,n) > max{A(k, m,,n), A(k, my, n)} 
and, in fact, it is easy to find numerical examples in which equality 
is false.t Thus it is easily verified that 


A(6,8,41)=1,  A(6,9,41)=4, A(6,72,41) = 5. 


We can, however, deduce from Lemma | that 


A(k, mms, n) < max{A(k, my, n), A(k, mo, n)} 
< 


and A(k,m,m,) < max{A(k, m,), A(k, mg)} \? 


provided that (m,,m,) = 1. 
We see that a knowledge of upper bounds for A(/,m) is useful if 
+ This difficulty appears to have been overlooked by I. Chowla in his paper 


‘On T'(k) in Waring’s problem and an analogous function’, Proc. Indian Acad. 
Sci., Sect. A, 5 (1937), 269-76. 
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we wish to calculate A(k), and we therefore study A(k,m) as fully 
as A(k,m). It does not seem unreasonable to suppose that 
A(k,m) = T(k,m), 

but we are not able to prove this. We can prove, however, that A(k, m) 
satisfies every inequality which ['(k, m) is at present known to satisfy. 

2.2. The numbers A’(k,m), 5,(k), and A,(k). We can simplify our 
problem by further considerations which apply somewhat similarly 
to A(k,m) and A(k,m). We call a solution of (1.2) primitive, if the 
greatest common divisor of m, x,, %9,..., %, is 1. We define A’(k, m,n) 
as the least s such that (1.2) has a primitive solution for suitable 
r = r(k,m,n), and we write 

A’(k,m) = max A’(k, m,n). 


Let p be any prime. We write 


k = p*k*, 
(k*,p) =1, 
(9+2 (p= 2) 


« = (p—1,k*), p—l= ed, k* = ek,. 
Where it is not clear to which prime we refer, we write 0,,, ¢,, €,, dy 
for 0, d, e, d. 

We define y,,(k) as the least value of s such that (1.2) has a primi- 
tive solution for m = p?, every n and r = s; and A,(k) as the least 
value of s such that (1.2) has a primitive solution for m = p?, every n 
and every r (0 <r <s). Similarly 5,(k) = A’ (k, p?). 

LemMA 2. If (m,,m.) = 1 and if there is a primitive solution of 
(1.2) for given n, 8, r when m = m, and when m = mg, then there is 
a primitive solution of (1.2) for the same n, s, r when m = my Mz. 

The proof is the same as that of Lemma 1. 

Lemma 3. If l > ¢ and x = y+zp!-?-|, then 

ak = ykt kzp!—9-lyk-1 (mod p’). 

This is Lemma 1 of P.N. 4. 

Lemma 4. If l > ¢, and if there is a primitive solution of (1.2) for 
given n, 8, r when m = p'-, then there is a primitive solution of (1.2) 
for the same n, 8, r when m = pi. 

The proof of this lemma is effectively the same as that of Lemma 5 


of P.N. 4 and so we omit it. 
3695.10 Oo 
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a4 > ol) — . 
Lema 5. A(k, p’) < A,(k) 


and A(k,m) < A'(k, ial per). 


pim 
The first part is obvious for ¢ < ¢ and it follows for t = ¢+1, 
¢+-2,... by successive applications of Lemma 4. If m = J] p*, the 
second part is obvious for t,, < ¢, (p|m); for larger t,, it follows by 
successive applications of Lemma 4, of Lemma 2, and of the trivial 
converse of the latter. 


Lemma 6. (i) Except when k = 4 and p = 2, 


A(k, p*) > X,(k). 


p 

(ii) Except when k = 4 and one of p,, po,... 18 2, 

A(k, pk pk ... p¥) > A'(k, por, poe, .++ piri). 

(iii) For all k and p, 

A(k, p*) > 6,,(k). 

Apart from the case k = 4, p = 2 we have always k>¢. If 
s = A(k, p*), then (1.2) is soluble for m = p*, all values of n and 
every r (0 <r <:s). This is true in particular for n = 1, 2, 3 
Since p* does not divide these values of n, the solutions must be 
primitive and so s > A,(k). 

Again let M = pk... p*, M’ = per, per. ... pe, and s = A(k, M). 
Then (1.2) is soluble for m = M, every value of n, and a suitable 
ry = r(k,M,n). If l<t< M’, one of t, t+’ is not divisible by 
pk, where 1<j <i, since p*/ M’. Hence either for n = ¢t or for 
n = t+M"’ the solution of (1.2) with m = M is primitive. This pro- 
vides a primitive solution of (1.2) for m= M’ and n =f, since 
M'|M. Hence s > A’(k, M’). 

The statement (iii) is a particular case of (ii) (with 7 = 1), except 
when k = 4and p = 2. But the fourth power residues (mod 2¢ = 16) 
are 0 and 1; hence 

5,(4) = 8, A(4, 16) = 8. 


Lemma 7. (i) Except when k = 4 and p = 2, 


A,(k) = max A(k, p'). 
t 


(ii) [fk A 4 orifk = 4and py, pag,... are all odd, then 


A’ (k, ph pi: ..+ phi) — max A(k, pp pe ... pt). 
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(iii) For all k and p 
5,,(k) = max A(k, p’). 
t 


(iv) A,(4) = 16, max A(4, #) = 15. 

The first three results follow at once from Lemmas 5 and 6. 

We see also that the congruence (1.2) has a primitive solution 
when k = 4, m = 2¢ = 16, and s = 15 (i) for all mn and 1l<r<l14, 
and (ii) for all nm 4 0 (mod 16) and 0 < r < 15, but not for r= 0 
or 15 and n=0 (mod16). It follows that A,(4) =16 and, by 
Lemma 4, that (1.2) has a solution with m = 24 and s = 15 except 
when r = 0 or 15 and nm = 16*l, where v > 0 and/ = 0 (mod 16). In 
the latter case there is a solution of 


l= 2t+...taf—at,,...—x#, (mod 2‘) 


for 0 < r < 15 and so a solution of (1.2), viz. 
l 


n = 16° = (2°x,)*+...+(2"a,)*—(2"a,,,)*—...—(2"a,5)* (mod 2%). 
Hence A(4, 2) < 15 and so 


max A(4, 2‘) = 15, 
t 


since A(4, 16) > A(4, 16, 15) = 15. 
Our calculation of A(k) is based on the following fundamental 
lemma which is a consequence of the above. 


Lemma 8. If we select a finite number of primes q, we have 
A(k) < max|A’(k, II @*), max A, (k)}. 
qd . pra 


Ifk #4, A(k) > max A’(k,m) > max3,(h), 


m p 

and A(4) > 8,(4) = 8. 

2.3. THEOREM 1. When k 1s odd, 

S2(k) = Ag(k) = yo(k) = 2; 
and, when k is even, 
8,(k) = 29+, A,(k) = 2°+8. 

When £ is odd, 6 = 0, d = 2, the kth-power residues to modulus 4 

are 0, 1, —1, and so 
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When & is even, Lemma 3 with p = 2,1 = 0+3, y = +1 gives 
(—1-+-42)* = (1+4z)* =1 (mod 20+2 
and we have (2z)k = 0 (mod 29+2), 
since k > 6+2. Hence 0 and 1 are the only kth-power residues to 
modulus 2°+2, and so 
A(k) = 2942, 
2.4. If —1 is a kth-power residue to modulus m, we have 


k 


ak — ak 
—H; = 


y; (modm) 


for some y;, and so, if (1.2) has a solution for any one r, it has a 
solution for all r. Hence 
A(k,m) = A(k,m) = I'(k,m). 
This is always the case if k is odd, since (—1)* = —1, and so we have 
THEOREM 2. If k is odd, 
A(k) = A(k) = T(f). 


Similarly, for any k, if —1 is a kth-power residue to modulus p?, 


we have 5,(k) = A,(k) = y,(h). 


Henceforth we suppose p > 2. Then there is a primitive root @ 
of p? and the kth-power residues to modulus p% are given by the 
d numbers Gute, G2ve Gav (2.41) 
(P.N. 8, pp. 521-2). Also 

—1= Gir*e-) — Gidr’« (mod p?). 
Therefore —1 is a kth-power residue, if and only if 
ap%e = 4dp% (moddp%e) 
is soluble, i.e. if d is even. Hence we have 
Lemma 9. Jf d is even, 
By(k) = Ay(k) = yp(h). 


If d is odd, the kth-power residues are again given by (2.41) and 
their negatives by 
Giap%(d+2/) (f= 1, By ccegt, 


i j ‘ 
that is, by Giv'+2), 
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where / runs through a complete set of residues to modulus d. The 
residue classes (2.41) and (2.42) together are given by the formula 

Gch, (2.43) 

where h runs through a complete set of residues to modulus 2d. 

Now let k = 2°’ and p—1 = 2», where k’ and 7 are odd. Since 
d and p are odd, o >7 >1. We write k, for 2"-*k’ and define @,, ¢,, 
e,, d, in terms of k, and p, just as we defined 0, ¢, «, d in terms of 
k and p. Then 6, = 0, ¢, = ¢, 

« = (p—1,k*) = 2°(n, 2-*k*), 
€, = (p—1, 27-°-1k*) = 27-1(2n, 2-°k*) = te, 
since 2-°k* is odd, and d, = 2d. Therefore (2.41) with k, in place 
of k shows that the k,th-power residues to modulus p® are 
Grah — Giv'ch (h = 1,2...., 2d). 

3ut these are exactly the residues (2.43). Hence the problem of 
representing a number to modulus p? by a sum of positive and 
negative kth powers is identical with the problem of representing 
a number to modulus p* by a sum of positive k,th powers. Com- 
bining this with Lemma 9, we have 


Lemma 10. If k = 2°’, p > 2, and p—1= 27», where k’ and » are 


odd, then ia [7 (oc < 7), 
4 Y,(27-°"k) (o >7). 

This reduces the calculation of 8,,() to that of y,(k) and the latter 
problem has been thoroughly investigated by Hardy and Littlewood. 

2.5. Upper bounds for X,(k). As we have seen, the primitive kth- 
power residues to modulus p? are the d numbers (2.41). Hence, if 
d =1, there is just one primitive kth-power residue, viz. 1, and 
if d = 2, there are just two, viz. 1 and —1. We deduce 

THEOREM 3. If p > 2 andd =1, 

A,(k) — Yp(k) a p?, 5,,(k) — 4(p?—1), 

while, if d = 2, 


A,(k) = yp(k) = 8,(k) = 1(p?—1). 


For larger values of d we shall prove the following: 


THEOREM 4. If p > 2, 0 = 0, and d >1, then 
A,(k) < max(e, 3) < k. 
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THEOREM 5. If p > 2, 0 > 0, and d > 2, then 
A,(k) < pe <k. 


Theorems 2 and 3 of P.N. 8 give the same results for y,(k). Hence, 
by Lemma 9, we need only prove Theorems 4 and 5 when d is odd, 
so that d > 3. Again the results for y,(k) cover the case r = s, and 
the case r = 0 is an immediate corollary. Hence, to complete the 
proofs of Theorems 4 and 5, it is enough to prove that, when m = p?, 
d is odd, and s = p%e, then (1.2) has a primitive solution for every r 
such that 1 < r < s—1 and for all n. In fact, if n 4 0 (mod p%), 
any solution must be primitive, since p? is a divisor of p*; while, if 
n = 0, the primitive solution 0 = 1*—1* is admissible, since 

l<r<s—l. 

Hence we see that Theorems 4 and 5 follow from the following 
two lemmas. 


Lemma 11. If p > 2, then 
A(k,p) << e<k. 
Lemma 12. If p > 2,d > 3, and 0 > 0, then 
A(k, p?) < p® <k. 


The proof of Lemma 16 of P.N. 8 can be generalized without 
difficulty to prove our Lemma 11, but it is simpler to use a theorem 
due to Davenport and I. Chowla,} viz.: 


Let {a} = {a,,d,...,a,} and {b} = {b,,b,,...,b,} be two different 
systems of residue classes to modulus h. No two a are congruent to 
modulus h and no two b are congruent to modulus h. Also (b;,h) = 1. 
Let | be the number of different residue classes c, to modulus h of the 
form a, or a;+-6;. Then 


_ 
Ll > min{h, m+n}. 


We take h = p. To prove that (1.2) is soluble with m = p, s = e, 
and a fixed r < s, we first apply the theorem r times in succession, 
starting with {a} = {0} and {b} = {x}, where x runs through a com- 
plete set of primitive kth-power residues to modulus p, and con- 
tinuing with the same set {b}, but taking as {a} the sum-set {c} 
obtained at the previous application. Next we apply the theorem 


+ Landau, Uber einige neuere Ergebnisse der additiven Zahlentheorie (Cam- 
bridge, 1937), Satz 114. 
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again e—r times, taking {b} = {—2x}, where x runs through a complete 
set of primitive kth-power residues to modulus p. We find finally that 


ed+1= p—1+l=p 
residues, i.e. all residues to modulus p, are representable in the form 


Since this is true for all r < e, Lemma 11 is proved. 


2.6. Proof of Lemma 12. We base the proof of Lemma 12 on a 
series of lemmas. 

We write N(s,r) for any » which can be expressed in the form 
(1.1) and we use a = No,r) 


(mod m) 


to denote that (1.2) is soluble for the particular m, n, s, r concerned. 
It is clear that 
N(81,7)+ N( (82,7 r,) = — N(8,+8:, r+. ), 

>and IN(s,r) = N(Is, Ir). 

In this section all congruences are to moduius p? except when 
it is stated otherwise. 

Lemma 13. If d > 3, then there is a kth-power residue c to modulus 
p* such that 

(p,c) = 1, ptli<c<(p—l1)p—l, c ~ ¥(p?+1). 

Since d > 3, there is at least one kth-power residue not divisible 
by p and not congruent to 1 or —1, say 2*. Then there is a number 
y such that zy = 1, so that a*y* =1 and y* #1. Hence, if z* =a 
and y* = b, where |a| < 4p? and |b| < }p*, we have |a| >1 and 
\b| >1. Also, if 

ja| < p—1, |b] < p—1, 

we have 1 < |ab| < (p—1)?, 
so that akyk = ab ~1, 
a contradiction. Hence either |a| > p or |b| > p. Let us suppose 
that |a| > p; then, if «* # }(p?+-1), we can takec =a= x*, where 
ptl<c<(p—l1)p—1. 

If 2* = 3(p*?+1), we have 

(a?)k = }(p*+2p?+1) = p*}(p?—1)+4(8p?+1) 
= 3p?+1). 
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Since d > 3, p > 5 and so 
p+1< }(3p?+1) < (p—I)p—1. 

Hence in this case we can take c = }(3p?+-1) # 3(p?+1). 

Lemma 14. If d > 3, then for any r (0 <r < p—1) we can find 
u = u(k,p,r) = 0 (mod p) such that 

up = N(p—l1,r) (mod >p?). (2.61) 

We write r’ = p—i—r. The case r = p—1 is Lemma 17 of P.N. 8, 
Hence we may suppose r < p—2 and?’ >1. We may further sup- 
pose r > r’ since otherwise we have only to change u into —u. Since 


» > 5, we have ea 
Pz . r>h(p—1) >2. 


Using the c of Lemma 13, we write 


p= Au pp =c+h. 
1 


Then (u—l)p+l<c<pp-l 


2<pxp-l, l<h< p—l. 
If h < r, we have 

pp = cth = N(A+1,h+1) = N(p—1,r) 
and pp + 0 (mod p); again, ifh > 7’ >1, ie. p—h <r, 

(l—p)p = —c+(p—h) = N(p—h+1, p—h) = N(p—1,r) 
and 1—yw # 0 (modp). The only remaining case is 
h=r’ =r= }(p—1). 

In this case 

(2u—1)p = 2c—1 = N(3, 2) = N(p—1,7), 


since r > 2 and?’ >1. This proves the lemma, unless 
2u—1=0(modp), ie. wp = $(p+l). 
But this implies c = pp—h = 3(p?+1) 
in contradiction to Lemma 13. 
Lemma 15. For any integers n and p (0 < p < p—1) 
np = N(e(p—1),ep) (mod p?). 
By Lemma 14 we can find wu + 0 (mod p) such that 


t., =. 
up = N(p—1,p) = Yyf— D> yf. 
i=1 i=p+1 








THE ‘EASIER’ WARING PROBLEM 201 
Since wu 4 0, we can find v such that wv = n (mod p). By Lemma 11, 
€ 
v= N(e,e) = DS x (modp). 
h=1 


Hence 


€ 
np = uvp = > xkup 
h=1 


€ p p-1 
=Ya(Yy¥— > yf) 
1 i=1 i=p+ 


h 


= Swi 


i=1h=1 i +1 
N(e(p—1), ep). 
We now deduce Lemma 12. We have to consider (1.2) with 
s = p®e, m = p*, and 0 <r < p%. We can write r in the form 
r= pote(pitp2P+psp”.-.+pep?), 
where 0< po Se, O0<p<p—l1 (¢=1,2....,6). 
For any n we have, by Lemma 11 and repeated applications of 


-1 


3 (x, yi)* 
h=1 


Lemma 15, 
= N(e,po)+p 
(€, pp) +N (e(p—1), €p})+-2 p* 
€, Py) + N(e(p—1), ep;) + pN(e(p—1), ep2)+n3 p® 


= N(€, py) +N(e(p—1), <p,)+pN(e(p—1), €p2)+ 
+-p?N(e(p—1), eps) +... +-p9-4N(e(p—1), epg)+1941 p9* 
: N(e+e(p—1)+pe(p—1)+...+p?e(p—1), 
Potepytepop+...tepgp?) (mod p?) 
: N(ep*,r) (mod p?). 
This is Lemma 12. 
2.7. We may sum up certain of our results for 5,,(k) and A,(k) as 
follows: 
Lemma 16. (i) Jf # >1, 
A,(2*) = A,(2°.3) = 2942, §,(2°) = 8,(2°.3) = 29+), 
(ii) If w is an odd prime and d > 0, then 
A,,(779(2—1)) = P41, 
§,,(7°(7—1)) = A,(409(2—1)) = 8,.(42°(7—1)) = 3(a® 1-1). 
(iii) In all other cases, 5,(k) < A,(k) < k. 
The results (i) and (ii) are immediate consequences of Theorems 
1 and 3; (iii) follows from Theorems 1, 4, 5, except when p > 2 and 
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1 or 2. If d=1 and k 4 p*(p—1), we have p%(p—1)|k and 
2p*(p—1). Hence 
Anh) = = pitt < 2(p?t1—p?) < 
If d = 2and k £ 4p*(p—1), we have }p%(p—1) | 
Hence A,(k) = (p9+1—1) < p(p—1) < k. 


We shall also require the following lemmas. 


k. 
kand k > p*(p—1). 


Lemma 17. Jf (m,,m,) = 1, m, > m, > 2 and m, is odd, then 
A’(k,m,mg) < max{}(m,+1), m,— 1}. 

Since 1 is a kth-power residue to modulus m, and to modulus mg, 
we see that (1.2) has a primitive solution for every n (i) with m = mg, 
s=m,—l, r= m,—2, (ii) with m = m,, s = m,—1, r=1, and 
(iii) with m = m,, 8s = 3(m,+1), andr = 3(m,—1) or r = 1 according 
to the value of x. In particular the residue 0 is always representable 
as 1‘—1*, Hence, by Lemma 2, (1.2) has a primitive solution with 
mM = M1 Ms, 8 = max{4(m,+1),m,—1}, and r = s—1 or r=1 ac- 
cording to the value of n. This is our lemma. 


Lemma 18. Jf 0’ > 0 and 6” > 0, then 
p? (p,—1) # pe “(pe—1), unless p, = po, 


ip? (p,—1) a p® (p2—1), unless p, = p. = 2, 
29° 3A p"(p—1), unless p= 3, 0 = 0” = 1 


27.3 4 4p" (p—1). 

The last two statements are obvious, while, if either of the first 
two were false, we should have p?’|p,—1 and p%’ | p,—1, so that 
Py < Pp, and p, < p,, a contradiction. 

In our numerical calculations we shall use the following lemma. 
We omit the proof, which is a simple adaptation of those of Lemmas 
11 and 12 of P.N. 8 

Lemma 19. Jf 6 = 0, s > 2 and 

p >. (e— 1)%8-Dis—2), 
then A,(k) < s. 

2.8. General results for A(k). We can now find the actual values 

of A(k) when k belongs to one of the following classes: 

(I)k=2° (@>}). 

(@}) k= ola 7 an odd prime, # > 0, and k& not in class I. 
) 


ps k = }(7—1), a an odd prime. and & not in class I or class II. 
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THEOREM 6.7} (i) If k belongs to class I, A(k) = 2°+, 
ii) If k belongs to class I1, A(k) = 4(7®!—1) > k+1. 
a ) If k does not belong to classes I or I1, A(k) < k. 
) If k belongs to class III, A(k) = 4(7—1) = k. 

Let k = 2° (# > 1); then 5,(k) = 2k by Theorem 1. If p > 2, then 

6é=0. If d =1, then 
(p—1)|k, p—l1l<k, A,(k) = p < k+1< 2k 

by Theorem 3. If d > 2, A,(k) <k by Theorem 4. Hence, by 


Lemma 8, 
A(k) < max{d, (k), max A p(k)} = 2k 


p>2 
and A(k) > 8,(k) = 2k. 
Hence we have (i). 

We say that k is exceptional with respect to 2 if k = 2°.3 (@ > 1) 
and that k is exceptional with respect to the odd prime z if k is not 
a power of 2 and if either 

k=7*(r—-1) (®@>0) or k= 4n*(z—-1) D1). 
Lemma 16 shows that A,(k) < k unless k is a power of 2 or k is 
exceptional with respect to p. We have to study A’(k, T] q*) for 
exceptional k, the product being taken over all primes q with respect 
to which k is exceptional. 

If k is exceptional with respect to one prime only, either 


k=2°.3 @>1), A(k, [I gs) = 8,(k) = 24 <b, 
o k= — (9 > 0), 8,(k) = (n9A—-1) < k, 
or k=4n*(z—-1) (@>0), §(k) = 4(w®A-1) > k+1. 
If k is wei with respect to each of two primes, either 
k=2?.3=n-1 (@>1), 
A'(k, 29+) < max{}(r+1, 2°41 ]} < k, 
k = 7*(x—-1) = w—1 (& > 0), 

A'(k, P®4H@w) < max{}(79®+1+1), o—1} < 
or k= 4n%(n—-1)=a—1 (8 >0), 8, (k) = 44-1) >k+1, 
and Aw(k) = wow = k+1, 


1 
2 
1 

2 


+ This includes two theorems stated without proof by Wright, Quart. J. of 
Math. (Oxford), 6 (1935), 266. 








204 W. H. J. FUCHS AND E. M. WRIGHT 


by Lemmas 16 and 18. Finally, by Lemma 18, the only / exceptional 
with respect to more than two primes is 6; and 
A’(6, 23. 32.7) <6 
by numerical calculation. Hence, by Lemma 8, 
A(k) < k, 

except when k = }7°(7—1) (® > 0), when 

A(k) = 8,(k) = 3(7®+1—1). 
To prove (iv) we have (iii) and 

A(M—1)) > 8,(4(7—1)) = (n—1). 

2.9. Calculation of A(k) for particular k. With the aid of the 
general results proved above, we have been able to calculate A(x) 
for values of k < 36 without excessive labour. Most of the remarks 
in P.N. 8 concerning the calculation of ['(/) apply, mutatis mutandis, 
to that of A(z). 

If k is odd, Theorem 2 enables us to obtain the value of A(k) from 
that of ['(k) calculated by Hardy and Littlewood. Theorem 6 gives 
the value of A(k) when k = 4, 8, 16, 32 (class I), 10 (class II), and 
6, 14, 18, 20, 26, 30, 36 (class III). The remaining values of k below 
36 are 12, 22, 24, 28, 34. 


re 


As an example of the calculation of A(k) for these k we give in 
detail that of A(12). We consider first for which primes q it is likely 
that ,(12) will be largest. By Theorem 4, we may expect to find 
these among the p for which 6 > 0 or for which « = k = 12, i.e. the 
primes /k+-1. Thus the least primes we have to consider are g = 2, 
3, 13. For the three primes 2, 3, 13, we have d = 1, so that 0 and 1 
are the only twelfth-power residues to modulus 16, 9, and 13. It is 


easily seen that 

A(12) > A’(12, 9.16) > A’(12, 9.16, 72) = 9. (2.91) 
Again every residue to modulus 9.13.16 has a primitive representa- 
tion with s = 9 andr suitably chosen among the values 0, 1, 3, 6, 8, 9. 
We note that either r > 6 or s—r > 6, so that there are always 
either six positive or six negative twelfth powers available. 

If p > 3 and p #1 (mod 12), we have 
A,(12) < max(e, 3) < 6 < 9 
by Theorem 4. Also A,(12) < 9 for 
p>11'4", i. p > 240 
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by Lemma 19. Hence, if we take the primes g of Lemma 8 to be 
2, 3, 13, 37, 61, 73, 97, 109, 157, 181, 193, 229, (2.92) 
i.e. the primes 2, 3 and those primes of the form 12/+-1 which are 
less than 240, we have 

maxA,(12) < 9. (2.93) 

pFa 
By calculation we find that y,(12) <6 for all the primes of 
(2.92) except 2, 3, 13; that is, (1.2) has a primitive solution for 


(i) m = TJ q, s = 6, r = 6, every n, and (ii) the same m, s = 6, 
a>13 
r = 0, every n. As we saw above, (1.2) has a primitive solution for 


m = 24.37.13, s = 9, every n, either r >6 or r < 3, and so, by 
Lemma 2, for m = [I q?, s = 9, every n, and suitable r. Hence 


: A’(12, TT a*) 9. (2.94) 
q 


< 
By (2.93), (2.94), and Lemma 8, A(12) < 9. Hence, by (2.91), 
A(12) = 9. 
Here we have an example of A(k) greater than max6,(k), for 
max 6,,(12) = 6,(12) = 8. . 


p 


TABLE 1 
k 2 ae 
A ; ; j 4 i¢ 


18 19 20 
18 20 


29 31 
26 15 29 5 64 33 


3. Upper and lower bounds for x(k). 


We write [ay,...)@]p-9 = 


to denote that 
j j 
>~a=— >a? C=1,3 
i=1 1 


i 


j j 

kml 2 k—1 
>a * + Sar-*. 
i=1 i=1 


It is obvious that (3.1) is equivalent to the identity (1.3) with 


j j 
C=G= i > at-1— > at}, 
i=1 i=1 











206 W. H. J. FUCHS AND E. M. WRIGHT 
Tarryt has remarked that (3.1) implies that 
[ay,...,4;,4,+h,..., +h], = [@j,-..,0j,0,+h,...,a;+h]},. (3.3) 
for every h; this provides an identity of the type (1.3) with k+1 
replacing k. There is a simple relation between the new C = C,,, 
and the C;, of the first identity. Thus 


Cray = (k+1) > {ak + (a;+h)*—aj*¥—(a;+h)*} 


= (k+ 1) tak 1_gq*k-}} 
1 


- (ene ® (3.4) 
In the special case 7 = k—1, we can express C;, in terms of the 
a and the a’ in a form which shortens its calculation. If 
A, = — } 4, A, = > 4,4 
are the symmetric functions of the a, and Aj, Ag,... the corresponding 
functions of the a’, and if 
j j 
S = > a, S; = ee 
t=] 
then Girard and Newton{ proved that 

S+4,45, it 2AS_.+...+1A, == 9 
and similarly for the A’ and the 8S’. Taking 7 = k—1, we see that 
A, = Aj, A, = Aj,..., Ap_o = Ay_s follow from (3.1). Hence 

= k{S,_,—Sy_1} = k(k—1)(Ay_y—Ay-4) 
k-1 k-1 ) 
= (—1)*e(k—1){ T] a,— [] aj). (3.5) 
i=1 v=1 
For k = 7 we start from [0,2], = [1,1],, use Tarry’s device four 
times taking h = 2, 3, 4, 5 in succession, and obtain 
[0, 3, 5, 11, 13, 16], = [1, 1, 8, 8, 15, 15],. (3.6) 
For k = 8 we use the solution of (3.1), viz. 

[0, 18, 27, 58, 64, 89, 101], = [1, 13, 38, 44, 75, 84, 102],, (3.7) 
discovered by Escott (and later independently by Phillips). For 
k = 9 we use the solution of (3.1), viz. 

[0, 4, 9, 23, 27, 41, 46, 50], = [1, 2, 11, 20, 30, 39,48, 49],, (3.8) 


+ L’intermédiaire des mathématiciens, 19 (1912), 200, 219-21; 20 (1913), 
68-70. 
t I. Newton, Opera omnia, ed. Horsley (London, 1779), 1, 182. 
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discovered by Tarry. In each of them j = k—1, and we can cal- 
culate C;, very simply from (3.5). 

From (3.8) we find solutions of (3.1) for k = 10 and for 12 <k < 20 
by using Tarry’s device with suitable h. For k = 11 we start from 
[0,3], = [1,2], and apply the device repeatedly. The values of h 
used are given in Table 2, as are the values of C which are found 
by the repeated use of (3.4). 

For k = 6 we can improve upon the upper bound obtained from 
the identities of the type (1.3) by using a slightly different type of 
identity, viz. 
l2abedfgx = (a°c+bdx)*+ (a°d—bex)*+ (b'c—adx)®+- (b°d+-acx)§— 

—(a°c—bdx)*— (a'd +bea)*— (b®c+-adx)*—(b°'d—acx)®, (3.9) 
where f = ct—d* and g = a**—}*4, This identity is due to Subba 
Rao. We have not been able to find similar suitable identities for 
higher values of k. 

If we take a, b, c, d positive integers, it is easy to see that (3.9) 


leads to v(k) < 2j+A(6, 12abedfyg), 
where 2j = 8 is the number of powers on the right-hand side of (3.9). 


TABLE 2+ 


Identity j 


(3.9) 4 12abedfg 
(3.6) 5 | 6.7.87.15? 
(3.7) 7 | 7.8.13.38.44.75.84.102 
(3.8) 8 | 23.3°.5*.72.11.13 
(3.8);h =9 0 | 9.10\C,| 
(0, 3], = [1, 2],;4 = 3, 4, &, 7, 8, i, 2 | 2%0.37.5*.73.11%.13.17 
13, 17 
Identity of k = 10; h = 17, 19 11.12.17.19|C,o| 
¥ » £=B;A = 6 6.13|C,.| 
= 13;iA= 4 4.14/C,5| 
14;h = 15 2 | 157|C;,,| 
15;h = 10 3 | 16.10|C,,| 
= 1l4;h 1, 7, 10 7.10.15.16.17|C,,| 
» k= Wk = 16 5 | 18.15|C,, | 
» & 18;h 23 56 | 19.23|C,,| 
b= 19;h = 16 | 6 20.16 /C\g | 





+ The identity for k = 6 was given by K. Subba Rao, J. of London Math. 
Soc. 13 (1938), 14-16; that for k = 8 by Escott, Quart. J. of Math. 41 (1910), 
152; that for k = 9 by Tarry, loc. cit.; those for k = 11, 17-20 by 8. Sastry 
and T. Rai, Indian Phys. Math. J. 7 (1936), 17-19; those for k = 12-16 by 
A. Moessner and W. Schulz, Math. Zeits. 41 (1936), 340-4. 
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It remains to calculate A(k, C;,). 
[k = 6.] It is easy to see that C, = 12abcdkl is always divisible 
by 7.8.9. Henze 


6 = A(6,7.8.9) < A(6,C,) < A(6) = 6. 
[k = 7.] We have C, = 6.7.8?.15. But since 
$= (p = 2), ¢ =} (p = 3, 5) 


for k = 7, it is sufficient to consider primitive solutions of con- 





gruences to modulus 2?.3.5.7. Since 
d=2 (p= 3), d=4 (p= 5), d@d=6 (p= 7}, 
8,(7) = 2, 
every number is a 7th-power residue to modulus 3, 5, and 7, and \ 


two powers are sufficient to represent any residue to modulus 


92 Qn 7 stricn 
2?.3.5.7. Hence A(7,C,) = 2. 


[k = 8.] Since 2¢ = 32 is a divisor of Cg, 
16 = 8,(8) < A(8, Cy) < A(8) = 16. 
[k = 9.] Since 27 = 3¢ is a divisor of Cy, 
13 = 8,(9) < A(9,C,) < A(9) = 13. 
[k = 10.] 5¢ = 25 divides C,, and 5,(10) = 12 since d = 2. Hence 
12 = 8,(10) < A(10,C,,) < A(10) = 12. 


[k=11.] If p|C, and p11, we have @=0, e=1, and 
therefore A,(11) <3 by Theorem 4. The 11th-power residues to 
modulus 11? are +1, +3, +9, +27, +81, and we find that 
§,,(11) = A,,(11) = 4. Also 11?|C,, and so 

4 = §,,(11) < A(11,C,,) < maxA,(11) = 4 
piC,y f 
by (2.12) and Lemma 7. 
Similar calculations determine the values of the remaining A(k, C;). 
In the third column of Table 3 we give the least M for which 


A(k, M) = A(k, C,). 
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TABLE 3 


Lower bound for | Upper bound for 
v(k) = A(k) | o(k) = 5+ A(E, ,) 
6 14 
7 14 
16 30 
13 29 
12 32 
11 
9 
6 
14 
15 
32 
6 
18 
4 
20 
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SOME GENERALIZATIONS OF A FORMULA OF 
RAMANUJAN 
By F. M. GOODSPEED (Cambridge) 
[Received 17 April 1939] 

1. IN a paper published in 1937 in The Quarterly Journal of Mathe- 
matics,* Hardy discussed the validity of two formulae of Ramanujan. 
These gave the cosine or sine transforms of a function defined by 
an infinite series, and a suggestion was made that they might be 
extended to the more general Watson transforms, with a formal 
indication of the method to be used. 

In the following, Theorem 1 is a more symmetrical and slightly 
more general form of the theorem for cosine transforms proved by 
Hardy. Theorem 2 is an extension to general Watson transforms, 
and Theorem 3 is a discussion of a special case closely related to 
a Hilbert transform. 

THEOREM 1. Let y(u) be an integral function satisfying the condition 

xlo+u) — = Ofedloi+Bitl), 
2tleiP(4|a|+ 42t) 


where B < 2, for all ao and t. Then 


S (—1)"x(n) an 

- 24"T'(4n+4) 

is convergent for 0 < x < e-4, and represents an analytic function L(x) 
regular for all positive x. Also 


vp 
| (: ) | ‘L(x)cos yx dx = My), 
7 
: 
=“. (—1)"yx(—n—1 
where My) = = : sine +) v 
7 2N+ 32 


for 0 < y < e~, and is the analytic continuation of this for all positive 


y. Also . 
9 
J6) | M(x)cos yx dx = L(y). 
7 
0 


* G. H. Hardy, ‘Ramanujan and the theory of Fourier transforms’: Quart. 
J. of Math. (Oxford), 8 (1937), 245-54. 
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‘In view of the similarity of this theorem to the theorem proved 

by Hardy, I will only sketch the proof. Suppose that 
—l<c<Q@, o<2< oe. 

Then an application of Cauchy’s theorem shows that 


ct+io 
1 ax(u)c“ du 
L(x) = —-— 
(*) Qari I; 24¢P(du-+-4)sin ww’ 


where x“ has its principal value. This is true at first for 0 << x < e-4, 
but since the integral is uniformly convergent in any interval 
0<8<2x<A< ©, it gives the continuation of L(x) for all posi- 
tive x. 

It follows that, if y is real and positive, 


9 in 
J6) [ L(w)eos ya dx 
7 
0 
_ |f2 i a ax(uja“ du 
PS J | cos Qi | aoe eared ” 


0 c—ia 


ct+io 1 co) 
1 mx(u) du 2 
= ——, [ x(w) : —| | 2“cos ya dx 
2nt J 24"D'(4u+-4)sin wun \z/ | 
m : 





c+io 


‘7 ax(ujy-“-! du 


Qari 2-Ku+DP(—Jw)sin zw 
c—t@ 
The inversion of the order of integration may be justified by a method 
exactly similar to that used in Hardy’s paper originally referred to, 
and so need not be dealt with here. 
Another application of Cauchy’s theorem gives 


y (—1) n ‘aiael 
J fu L(x)cos yx dx = — : S1n- aM ae 


=> 1)"x(—2—}) ym — My). 


Tn +4) 
This again holds in the first place for 0 < y < e~4, but since the 
integral representation is uniformly convergent for 
0<8<y<A<o, 


the final result holds for all positive y. 
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In the above, x(u) satisfies conditions symmetrical about the origin, 
so we may replace y(u) by x(—u—1). This merely reverses the roles 
of L(x) and M(x), and so we obtain the reciprocal formula 


9 ve) 
J6 | M(x)cos yx dx = L(y). 


2. Examples. I select two examples which are not covered by 
Hardy’s theorem. Let 
2T(Ju-+4) 

T(n+1) 
The conditions of Theorem 1 are satisfied with A = 0, B = 31, and 
we obtain as reciprocal functions 


L(x) =e, M(x) = /Oae 


In Hardy’s theorem, as in Ramanujan’s work, both L(x) and M(x) 
are analytic for all real x and so, being even, are expressible as power 
series in x?. This restriction involves an additional condition on y(u). 
Here L(x), if defined by evenness for negative x, is not analytic at 
the origin. Hence the above theorem is more general than Hardy’s 


x(u) = 


result. 
Hardy states his conditions in terms of ¢(w), where 


x(u) o(u) 


2ieP(du+4) P(w+1)’ 


and ¢(w) = 0 for positive odd integral u. Apart from this restriction, 
his result is equivalent to Theorem 1. 
Another example is given by x(u) = 1. In this case we obtain the 


self-reciprocal function 


x 


L() = M@) => (—2P ns 
“, "Pn +4)” 
which has been given by Hardy and Titchmarsh.* In the paper 
referred to, this function is given in the form 
oo 0 
1—2ae“ | e- dt = [ e-em dw 
x 0 
which on expansion in powers of x yields the series obtained here. 
* G. H. Hardy and E. C. Titchmarsh, ‘Self-reciprocal functions’: Quart. 
J. of Math. (Oxford), 1 (1930), 210, example 4. 
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3. THEoREM 2. Let p(u) be the reciprocal of an integral function 
taking real values on the real axis. Let x(u) be an integral function 
satisfying the conditions 


x(o+it) = O(e4otBitl) x(—o+it) _ = O(eCo+Ditl 
p(o+1+it) "  p(o+it) 
where B<a, D<7, fora >0. Then 
be) = 20) 
p(i—s) 
satisfies the equations k(s)k(1—s) = 1, and |k(4+-it)| 


1. If 


a $+io 
K(x)_ 1 f R(s)a* 5, 
= 2m l—s 
t—iw 

where the integral exists in the mean-square sense at infinity and K,/x 
is L*(0,00), then és 
ya 
converges for 0 < 2 < e-4, and represents a function L(x) regular for 
all positive x and of integrable square in (0,00); and 


£ (SOY 12) dx = My) 
dyJ «2 
0 


for y > 0, where 


— (—1)"x(—2—1) 

ee gs - way y 

for 0 < y < e-©, and is the analytic continuation of the power-series, 
for all positive y. Also 

. | KY) y(n) da = L(y) 


dy , 
0 
for y>0. Thus L(x) and M(x) are Watson transforms with kernel 


K, |x. 


x 


We have, by Cauchy’s theorem, 


ct+io 
— (—1)"x a 1 ax(uja” du 
lee em (n+1) “ pio! }  p(w+1)sin wu 


c—tx 
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if —1<c< 0. As usual, the integral representation holds at first 
for 0 < x < e~4, and by analytic continuation for all positive x. It 
is plain that L(x) = O(a*) both for large and for small x, and hence 
taking in turn c > } and c < }, we see that L(x) is L*(0, 00). 

It follows that 


L(x) dx = —— : 
: (x) 9 J  sinzup(u+1) « 


© 


A t+i 4 , 
Aa(ey) : a at a Hey) dx 


. 


r) 4—ix 
for y > 0. But, from the theory of the Mellin transform, 
A 


a 


> 


Cy 


as A oo and 6+ 0, converges in mean square to —y~“k(u-+1)/u, 


and also 
x(u) 


sin zu p(u+ 1) 
is L?(—4+-i00, —}+7i0). Hence, by the usual theory of mean con- 


vergence, 


B ay(ujy-"k(u+ 1) du 
wsin zu p(u-+ 1) 


% —Uu aT = n-1 bee 
™x(w)Yy du = (—1)""x( 2) yn, 
J usin ru p(—u) i np(n) 


The final series converges for 0 < y < e-©, and as usual represents 
a function regular for all positive y. Hence 
8 

af Ki (xg 

—— | L(x) i( YD ay = M(y) 

dy J x 

0 

for y > 0. The reciprocal relation can be obtained at once from the 
general theory of Watson transforms or by replacing x(u) by 
.(—u—1) in the relation already obtained. Also 
x s 

<8 x 

{L(x)}? dx = | {M(x)}* dx. 


0 0 
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4. Examples. (i) Let p(u) =1. Then 


=F (—1x(n)e 


and K, = 0 for x < 0, K, = 1 for x >1, so that 


i = dz = 17) = 2 (—1)"x(—n—1)y". 
ily 

The series 2 (—1)"x(n)y-"-1 converges for large y, and the series 
> (—1)"x(—n—l)y" converges for small y. These series may not 
converge i any common values of y, but represent the same analytic 
function. Suppose, for example, that x(w) =1. Then 

i. 2 

~ +e’ 


and this function is expansible as } (—1)"2-"-1 for x >1, and as 
> (—1)"a" for x <1. 
(ii) Taking p(u) = 2K“-DP(42), 


/ a) (—1)"x(n) ‘ 
L(x) = > eT nth) 
- > 2 


and Sf 14) ax = My). 
dy x 
0 
This is the integrated form of Theorem 1. 
(iii) If p(w) = cosec 4zru, 


= ¥{ (—1) n x(2n)a?”, . —2n—1)a2", 
0 
and : s [ Le) lo 


andy) 2x 
0 


Formal differentiation under the integral sign yields 


o 


> f 2 dx = My). 


Cg —x*y* 
0 


I do not attempt to justify this differentiation directly, but show 
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that the result is true by an application of the method used in the 
proof of Theorem 1. 
5. THEOREM 3. Let x(u) be an integral function satisfying 
x(o-+it) = O(edlei+Bitl 
with B < 43n, for all o and t. Then the series 


2) 


y3 (—1)"x(2n)a2” 


converges for e~4 < x < e4, and represents a function L(x) regular for 


all real x; and = 


Ra lx = My), 


@ 
where = > (—1)"x(—2n—1)y” 
0 


for —e-4 < y < e~4, and is regular for all real y. Also, 
r b x 
[ = dx = L(y). 
1—z*y* 
0 
As in Theorem 1, we have 
L(x) = ¥ (—1)"x(2n)a2" 
0 
c+io 
y nl aul “U 
7x(u)sin 37(w-+ 1)a ie 
sin 7u 





-4, and the 


for —1 <c < 0, the series converging for —e-4 < x < e 
<A<o. Suppose 


integral converging uniformly for 0<é8<2z<A 
now that y is fixed and y > 0. Then 
y*- A 
i a 
y*+n 
yo 


ct+io 


' “Ss 7 

1 [ mx(u)sin $7(u+1) du [ [ x“ du 
a6 . Tv 24/2 
270 J sin 7u j ; 1—wx*y 


c—iw 5 yt’ 


for y > 0. Also, for small 2, 


et 
= = Q(z") 
1—z*y* 


yu 
* 5 = O(x°-2) 


and, for large x, 3-5 
1—z*y? 
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uniformly throughout the range of integration of u; and so it 
follows that we may replace 8 by 0, and A by o0 in the above 
relation, so that 


(/ + {es 


y "+n 





a7x(u)sin $n(u+1) du ae dx 
sin 7u f+ [ 1—2y? 


y7+7 


1+y9 
9 a” dx | ~u-a| t” dt 
1—z2y?| y 1—#’’ 
7% 1-vq 
. t¥ 1 l {u—] 
and, since a eee w= 
7" 1—# ~ 20—T 2H 
+yn 1+ 
[ "yu dt 1 dt 
7 ifs 1++# 
1—yy 1—-y 1—-yy 
Hence 


yt+n 
x” dx 
| i < ante [I ~ t/ dt < On-+Dnlul, 
lgt—@ 1-yy 


where A, B, C, D are independent of the imaginary part of wu. 
Therefore 


ct+io y*+7 

[ x(uw)sin 42(u+1) du x” dx 
j sin 7u 1—x*y? 
c—in y*-9 





is convergent and tends to zero as 7 > 0; and 
oo ctio oo 
L(x) dx _ 1 | 7x(u)sin 37(w-+1) du 2 =f x” dx 
— -s 


sin 7u l1—a*y? 


277 
c—ia 


ct+ia 
1 mx(u)y—“—* du = 
th aoe yi 
2nt J 2sin}a(u+1) Y) 


cC—12 


for 0 < y < eA, 
As before, since the condition on x(w) is symmetrical about the 
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origin, we may replace x(u) by x(—uw—1), and thereby obtain the 
reciprocal formula 2 


2 ( M(x) dx 
= = L(y). 
7 | 1—2*y? ) 


Example. Let x(u) =1. Then 
] 


1-2?’ 


and we obtain a function self-reciprocal with respect to the kernel 
> 3 
na 1—2?" 
Generally, any function self-reciprocal with respect to the kernel 
K,/x, and obtainable from Theorem 3, will be of the form 
c+ia 


1 f wx(ujav du = s ( 


a oe =—— > (—1)"x(2n)a?” 
2m 2 sin $aru a x ‘ 
c—io 
where x(w) is an even function of u+4. For instance, if 
x(u) = (w+3)*", 


we obtain the sequence of self-reciprocal functions 


at om a | 
dx ' 2) 1+22° 








HYPERGEOMETRIC PARTIAL DIFFERENTIAL 
EQUATIONS (III) 


By T. W. CHAUNDY (Ozford) 
[Received 21 March 1939] 


1. The equation y"f(8)V = x™g(8’)V 
IN two previous papers under this title (2, 3) I have discussed the 
solution of those partial differential equations of the form 


io oS 
fO\F OW =ayg) G0 (8,3'=22,y5) a) 


in which the operators of ‘dominant’ order are found on the same 
side of the equation: that is to say, if p, g, P, Q are the degrees of 
i. 9, F, G respectively, then 


pP>4 P>Q, (2) 


(or, of course, both inequalities might be reversed). I now consider 
certain equations of the same form (1) under the opposite conditions 
p>q@ P<Q. (3) 

I had previously proposed to distinguish these two classes of equa- 
tion as hyperbolic and parabolic, but I have since realized that this 
would be a misleading extension of accepted terms. The real dis- 
tinction (as I see it) arises in this way. Equations of the first type 
have for their characteristics the two families of axi-parallels 
x = constant, y = constant, i.e. the characteristics are always equal 
in number to the independent variables. Equations of the second 
type fall into two classes. If the aggregate orders of the operators 
are the same on both sides of the equation, i.e. if p+P =q+@Q 
in (3), then the equation will have other characteristics instead of 
(or in addition to) the axi-parallels and so, in general, will have 
more characteristics than there are independent variables: the 
characteristics will be ‘excessive’ as I have elsewhere proposed to 


say. Thus the equation 
s foes . — eV 
8(8—1)(8—2)V = ay-38’(8’—1)(8’—2)V, ie. in aye 
has the three families of characteristics 


x+wy = constant (w* = 1). 
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On the other hand, in the equation 


» =. ar - oF 
—1)V = a2y-18') ig ee ee 
5(6—1)) xy1d8'V, i.e @a% = By’ 


where the orders differ on the two sides, the characteristics are the 
single set of axi-parallels y = constant, i.e. the characteristics are 
‘deficient’, and the equation is parabolic in the proper sense of the 
term. 

In this paper I shall still exclude the properly parabolic case and 
consider only the equations in which p+Q = P-+gq, when (as I have 
said) the characteristics are, in general, excessive. An equation 
evidently of this type is* 

FE)V = xry’g(8')\V (4) 
in which f, g are polynomials of equal degree m. By a suitable trans- 
formation x’ = x", y’ = y* we can reduce it to the simpler form 

y™fO)V = am g(d')V, (5) 
in which we may further suppose that the leading coefficients in f, g 
are unity. In this case the terms of highest order will be 
y™ 8(5—1)...(8—m+1)V —a™ 8’(5’—1)...(8’—m+1)V, 
i. my amv mV 
exam oy” , 


and the characteristics are therefore the m pencils of straight lines 
x+wy = constant (w” = 1). 

I shall use a technique similar to that with which I have previously 

discussed a class of equations with ‘excessive’ characteristics.f 


2. The solution of Cauchy’s problem 

As usual we suppose given an arbitrary point P with coordinates 
(X,Y) and an arbitrary base-curve. The solution of ‘Cauchy’s 
problem’ consists in determining the value V(X,Y) at P of any 
solution V(«,y) of the given equation in terms of its value and that 
of certain of its derivatives along the given base-curve. 

Draw through P the pencil of the m characteristics 

x—X+w,(y—Y)=0 (r=1.,..., m) (6) 

* This, of course, does not represent the most general equation of the type. 
The general equation, with operators in 6 and 8’ on both sides, will have 
axi-parallel characteristics in addition to those of the type «+wy = constant. 


+ See (4). 
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cutting the base-curve in the m points A,,..., A,,, and take any con- 
venient origin O on the base-curve and any convenient curve OP 
connecting O to P. 

I then postulate a set of m ‘fundamental solutions’ (or ‘Green’s 
7 U(X,¥;a,y) (r=1...., m) (7) 
each associated with a characteristic PA, and having the following 
properties: 

(A) they are solutions, in variables (X,Y), of the given equation 
(5) and, in variables (x, y), of its adjoint equation; 

(B) they satisfy the ‘null condition’ 


m 
> U,=9; 
r=1 

(C) along the characteristic PA, the associated solution U, vanishes 
with all its derivatives in x,y of order less than m—2: its derivatives 
of order m—2 have the values 

g~-87/, a witl x a y "a yXY)- _a-D 
da™—2-aoyt sm \X) \Y 
along this characteristic, ma, mb being the coefficients of the terms of 
order m—1 in f(8), 9(8’). 

Observe in (C) that, since every derivative of U, of order m—3 
vanishes along PA,, we necessarily have, on differentiating along 
this line, gm- 207, gm- 207, 
Oram 2—19y2 exm- 3—doyatt 

m—277 
sera a wt 6, (8) 
for some 6,; the second part of the condition (C) serves merely to 
fix @,, defining it as 


m 


= 0, 


which gives 


I now make the convention that A, \’ denote the operators 6, 5’ 
acting only on V and its derivatives, and —y, —p’ the same opera- 
tors acting only on the U, and their derivatives, so that, on mixed 
operands in U,, V, the full operators are 

5 = A—p, o’ = N’—yp’. 
Thus, if V is any solution of the given equation (5), 
[y"fr)—a g(A’)]U, V = 0; 
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and, since, by (A), every U, is a solution of the adjoint equation 
ly f( = —a g( —$)|U = —_ 0, 
we have also Ly" f(u)—a™ g(u’)|U,V = 0. (11) 
By (10), (11) 
| Ayn MO) — yy’ ) gmd A) —g9(u! ‘l0.V = 0 
i.e. 0) ~- 3 yn) gy] 3 | )—9v My. v| fan 0. 
‘ A—p N’—p 
Thus over any closed area 


| ea"! . 0, 
xy 


and so, by the Green-Stokes theorem (in 6, 8’ as fundamental 
operators), 

| ( gd g(r’ Ie uy dx - yf) oy V dy)\ = 0 

fa F =p y| 
round the boundary of this area, i.e. round any closed curve. 

Take the closed curve to be the curve OA, PO made up of the are 
OA, of the base-curve, the segment A, P of the characteristic corre- 
sponding to U,, and the.arbitrary arc PO; and sum in r. Then, along 
the are PO, which is common to all the m integrals, U, appears as 


> U,, which vanishes in virtue of the null condition (B). Thus the 


part of the summed integral corresponding to the arbitrary curve 
PO vanishes, and we have left 
m P A, , ’ . 
= [+ [ {(ame@—a) oy iy ff) 7 y\\ Y} — 9 
= "ff oe ee 2 i a MR 
4, O 
(12) 

The second integral in the summation is taken always along the 
base-curve, i.e. in terms of Cauchy’s problem, it involves only 

y8 ] : 
‘known’ values of the unknown V and its derivatives (apart from 
the functions U, 


reduction of the first integral 


supposed known). Thus we need consider only the 


P 
[ (= g(A )— 9H) U, V dx + ra U, v) dy) (13) 
J Np ax A—p yJ 


taken along the characteristic A, P. 
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3. Reduction of the integral along the characteristic 

In considering the reduction of (13) I shall, for brevity, omit the 
suffix r from the associated A,, U,, w,, 6,, restoring it when I come 
to substitute the reduced value in the summation (12). Now, in the 
first place, in view of (C), we need retain along this characteristic 
no derivatives of U of order lower than m—2, i.e. we retain only 
powers of pu, »’ of orders m—1, m—2. To this order we have 


a UV = (w+ p+ ma p"™)UV 
.— 
— V[(—8)"-1+ ma(—8)"-2]U + (8V)(—8)"2U, 
and so for g. Again, sufficiently to this order, 
§”-1 — §(8—1)...(8—m-+ 2)+4(m—1)(m—2)8(8—1)...(6—m-+3) 
am—1 am—2 


— »m-1 jt 1(m— _ 9\ym—2 
=a Perna fe $(m—1)(m—2)x = 
and so for 8’. Thus, along PA, the integrand of (13) is effectively 
(after rearrangement) 


(- ay{V | er-1U 4(m—1)(m—2)—ma e"™"U 


dy T dy™-1 dx - . eym-2 dy = 
_ 3(m—1)(m—2)—mb aw eV mw _W om8U | 
' y oy*-* 


oxm- 1 x 


— ——, dz}. 
ox eam-2 y oy dym-2 | 
(14) 
Now differentiation along PA is effected by the operator wé/éx—é/éy. 
Thus differentiating w”-10 along PA we get, by (8), 
am—1]T7 am—1{7 2 a 
wma ¢ U —wm"-a-1 ¢ U o-—_ wml wl? 5 
da™—I-loyd da™—-I-2G +1 Cx oy 
Summing this in g from 0 to m—2 we have, since w” = 1, 
ty gn~tiy 06 06 
7 —o-; = (m—1)|——w* — }. (15) 
eam-1 oy™ -1 ex ey 
Then, in (14), taking the first pair of terms in the coefficient of V 
we have 
er -1UV om-1[y 
v Aqgym-l dx 
dy” 
(Fa or -1U7 


oam-1 hit oy™-1 


Caml 
)ay, since dx+wdy = 0 along AP, 


26 
cy 


= m— —w 
Ox 


Jay, by (15). 
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The second pair of terms in the coefficient of V can be written 


aR a a gg 
le al! ast lady 





x wy 
00 00 6 86 
on ee eA nin ee 
m (" w =) dy —(m (2 =| dy, 
since, by (9), 
28-188 _ fa—}(m—1) | b—3(m—1)\p 
Ox oy \ x wy | 


Thus the coefficient of V in (14) reduces to 
-(F-« 3F) ay —(m—1)(2— =) dy. 
6a dy wy 
The whole expression (14) can therefore be written 


__ gs\m—-1 
(—zy) |—¥ (0-2) ay — (m—1)V (¢- ;)ody — (wo —F) au} 
2-2 dx by 


Ww | Ox Oy 


OR ee =("(e dx +5 dy) + (m—1)0V (S+)+ 


WwW 


ve a 


—])m-1 
aan ( 1) d(x™—ly"—10V). 
wW 


The integrand (14) is thus an exact derivative along PA, and so 
the integral (13) reduces (with restored suffixes) to 


__})m-1 
( 1) [am—ly 1,;m— 19. vy 


WwW 4, 


me (—1)"-1 x a+4(m—1) y b+4(m—1) a P 
ra m = Y v| » by (9), 


A, 


- (—ip-*} 2 a a+k(m—1) y b+4(m—1) ~ 
foeen ly} 


Hence finally from (12) we get 
V(X,Y) 


— yn Sf ffm =I) 7 7] 8 [mf Fl) dy) 
=a 5 ffl: =, uv |= +|y 7 u,v |) + 


= a \a+h(m—)) y b+4(m—1) 4 
mal D3 Vi - (16 
+ m Zs (=) (¥) | wit 


nM 


r=15 
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The terms to be evaluated at A, can be replaced, if it is preferred, by 


A, 
1 me ae \a+hm—1) y b+Hm—-1) ( (QV a+43(m—1) ! 
> / G yet set 
r=10 


ae. b+3i(m—1),, ap \A+Hm—D) / y\b+Hm—1) Oo 

Lf = Vid = wo V > 
so that the whole expression on the right of (16) is then composed 
of a sum of integrals from O to A, together with a term at O. 

In either form (16) provides the solution of Cauchy’s problem,* 
for it gives the value of any solution V at the arbitrary point (X, Y) 
in terms of the values of V and its derivatives along the given base- 
curve, so long as we know the set of fundamental solutions 

U(X, Y;2, y). 
To the problem of determining this fundamental set we now turn. 
4. The fundamental solutions 


To avoid excess of algebraic detail I take the third-order equation 
as sufficiently illustrative. As a further simplification I note that 


the substitution V = aty-*y’ (17 4) 


transforms the given equation (5) into the form 
fy" f(8—h)—a” g(8’—k)}V’ = 0, 
of which the adjoint is 
{y™ f(—8—h)—a™ g(—8’—k)}U’ = 0. 

Hence the substitution 

U, = ff XT, (178) 
simultaneously applied to a set of fundamental solutions of (5) pre- 
serves the property (A) of §2. The null condition is evidently in- 
variant for such a substitution, and, lastly, the terms of order m—1 
in f(8—h), g(8’—k) have coefficients m(a—h), m(b—k), and so the 


* Or, more strictly, a possible solution of Cauchy’s problem. For, if m > 2, 
not all the characteristics are real. Thus some A, will be complex and the 
integrals will be, fundamentally, integrals in a complex plane, their values 
possibly depending on the path of integration. But the argument appears to 
show that (16) remains a solution for all such values. The underlying problem 
is then how far in such a case Cauchy’s problem is ‘correctly set’, ie. how 
far the data are adequate to define a solution. That is a problem I do not 
attempt to discuss. 

3695.10 Q 
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substitution also preserves the property (C). Thus a set of funda- 
mental solutions remains a set of fundamental solutions. 
I use the substitution to remove the absolute terms in f, g, so that 
the differential equation may be considered in the form 
y3(8+3p,+1)(8+3p.+2)V = 238'(8’+39,+1)(8’+3q.+2)V, (18) 
and I begin by considering the case in which p,, ., 91, do are all 
positive integers.* The adjoint of (18) is 
y® 5(8—3p,—1)(6—3p,—2)U = 238’(8’—3q,—1)(8’—3q,—2)U. (19) 


To obtain its solution I introduce the operators 


Ai(Py) = (6—1)(6—4)...(6—3p,+2) 

A?(po) = (6—2)(6—5)...(6—3p,+1) (20) 
Ai(q) = (8’—1)(6’—4)...(8’ —3q, +2) 

A2(q,) = (8’—2)(8’—5)...(8’—3q,-+1) 


the factors in each product descending by differences of 3, and 
substitute in (19) 
U = AX(pyAP2)AG(G)AG (Qs) W. (21) 
The equation can then be rewritten 
(8—4)...(8—3p,—1).(8—5)...(8—3p.—2).(8’—4)...(8’—3q,—1) x 
x (8’—5)...(3’—3q.—2)[ y2 8(8—1)(8— 2) —a* 8’(8’—1)(8’—2) ] W = 0 
and is therefore satisfied, if W is any solution of the differential 


equation 
[y®3(8—1)(8—2)—a8’(3’—1)(8’—2) ]W = 0, 
SW Ww 
i.e. of the equation wd : ae y = 0, 
dx oy 
Thus, from (21), a possible solution of (19) is 
U = Ah(py)AZ(P2) Ai(G)AF(G2) $(@+ wy), (22) 
where w* = 1 and the function ¢ is arbitrary. 
We now need to distinguish the signs of p,—p,. and q,—q». It is 
sufficiently illustrative if we suppose 
P2 > Pr h > +l. (23) 
Substituting Vi = o-BPat2y—Baty’ (24) 


* The partial differential equation has then an analogy with the ordinary 
differential equation discussed in (1). 
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in (18) we can write that equation as 
y® 3{8—3(p.—p,)—1}{8—3p,—2}V’ 
= #8 8{8—3q,—1}{5—3(q,—g.—1)—2}V’. 
Since po-—P,, Y¥;—2—1 are positive, as well as p,, q;, in virtue of the 
assumption (23), this equation is similar to (19) with changed para- 
meters. Thus, by (22), it has a solution of the form 


V’ = Al(p,—p,)A2(p2) AMq,)AX(q—G—1) Hw+ey). (25) 
We transform this to variables X, Y in accordance with the terms 
of (A),* and then combine (22), (24), (25) in the form 
U,, = $X—-BP2+2)(wV)-Oa+VAN (ps AZ (po) AF (G1) AF (Go) X 
w hte 2(y \Al A2/ {(x—X)+o(y—Y) osu 
x Ak(pe—Py)AX(P2) AY(G)AF-(G1—G2— 1) (3p,+-39,)! . 
(26) 
This, then, is a solution, in variables x, y, of (19), and, in variables 
X, Y, of (18), and so condition (A) is satisfied. 


5. Behaviour along a characteristic: the null condition 
Now along the characteristic 

(x—X)+w(y—Y) = 0 (27) 
the operand {(~a—X)+w(y—Y)}§?:+8% and its derivatives of order 
less than 3p.+3q, all vanish. The order of any operator A(h) is given 
by its parameter h, and so the eight operators in (26) have an aggre- 
gate order 3p,+3g,—1. Thus U,, vanishes along its associated 
characteristic (27). 

In 8U,, the aggregate order of the operators equals the degree of 
the operand. Thus, along the associated characteristic, 5U,, has a 
value which is not zero and which can be found by retaining only 
the operators of highest order. This gives along the characteristic} 

8U,, = 4X ~BP2+2)(qy V )—-BaitDyPit Path (qyy )d1+d2 X 2P2t+Pi(qyV )*di-de—1, 


i.e. oU,, a= dw Prt Pry td2X Pit P2t DV Gi tdet2), (28) 


Ox 


* The transformation requires us to define Ay, Ay with such a convention 
of sign that 
Ak(p) = (—8x+1)(—8x+4)...(—8x + 3p—2) 
= (—)?(8x—1)(6x—4)..-(8x—3p+2), 
and so on. 
+ We must not overlook the conventions of sign inherent in Ay, Ay. 








228 T. W. CHAUNDY 


which is the form required by (C), since, in (18), 


a= p,+po+l, 6b=q,+¢@+1, 4(m—1) = 1; 


so for dU,,/éy. 

To prove the null condition (B) we observe that U,, can be ex- 
panded in the form A+ Bw+Cw*. Since B, C disappear in the 
summation > U, we have to show that the term in U,, independent 
of w is absent, i.e., from (26), that the coefficient of w is zero in 


AH PyAZ(P2) Aj (G)AG(de) X 

x Ak(P2—Py)AK(P2) AV (GAH —d2— 1N(#@—X) +-o(y— ¥) rea, 
Now the coefficient of w in the expansion of {(w—X)-+-w(y—Y)}5?2+8u 
is composed of terms of the form 

(a—X)%*-1(y—Y)%*+1 (h+k = pot+q,) 
(of course, with appropriate coefficients, which I now ignore). The 
terms in the expansion of (z—X)*"-1 have one of the three forms 
(i) 34 X3b-1 (ii) 34-1 X 3b (iii) 34-2 X 3b+1 

where a+b=h. Of these (i) is annihilated by A%(p,) unless 
b> pot; (ii) by A2(p.) unless a > p,+1; and (iii) by A2(p,) 
unless a > p,+1, and by Ak(p,—p,) unless b > p,—p,. ‘Thus no 
term survives in (e—X)?"-1 unless, in every case, a+b > p,+1, 
i.e. unless h > p,+1. 

Again, the expansion of (y—Y)***+! consists of terms of one of the 

three forms 

(iv) yey ri (v) getty, (vi) gts e-1 
where a+b = k. Of these (iv) is annihilated by A;(q,) unless b > q,; 
(v) by A,(q,) unless a > q,; and (vi) by A,(q) unless a > q,, and by 
Ay(¢:—2—1) unless b > g,—qa. Thus no term survives in (y—Y)**+4 
unless, in every case, a+b > q,, i.e. unless k > q,. 

Hence one or other binomial in the product (~—X)*"-'(y—Y)3*+1 
is annihilated by the product of the A operators unless simultaneously 
h >p.+1 and k >q,. But this is impossible since h+-k = p.+q,. 
Thus the functions (26) satisfy the null condition (B) as well as the 
other conditions (A), (C) predicated of a set of fundamental solutions. 

This completes the solution of Cauchy’s problem for the third- 
order equation (18) in the case when p,, ps, 9g), do are positive 
integers.* 


* Or, sufficiently, when they are integers of either sign: a transformation 
of the type (17) reduces the more general case to the case we have considered. 
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This is the case in which the fundamental solutions are poly- 
nomials. In preparation for the more general case in which they are 
transcendental functions I put the expression (26) into the form of 
a contour integral. 


6. The fundamental solutions as contour integrals 
Since {(a—X)-+w(y—Y)}*?2+51/(3p.+39,)! is the coefficient of 
{972+3% in the expansion of exp{(a—X)t-+-w(y—Y)#}, we can write the 
fundamental solution (26) as the contour integral 
(0+) 
pp XOrm(w)-o0rem | AL pA IAHa)AH a) 


6772 
x Ak (po—Py)AR( pa) AY(G)A-(G1—G—1) X 


x exp{(~a—X)t+w(y—Y)} os 


{3P2+3qi+1? 


which I rearrange as 
(0+) 
Us = gj [ AMPDAR Met. Aq) AHae™ 
x (Xt)-Gr2+PAL (Po—p,)A 2(p Je -Xty 
x (w¥t)-OuPAL(q,)AP(G—d—le@F!.# dt. 


Let us write this 
(0+) 


U,, = 5 | £(at)n(wyt)E(Xt)H(wY te dt, (29) 


where E(at) = AMp,)A%p.)e% 
E(Xt) = (Xt)-GP2+?Ah(p.—p,)AX(pa)e*",_ ete. 
Thus &(2) = Al(p,)A2(p,)e”, since the operators A, involve only 4; 


and so for the other 6 operators. 
Consider the ordinary differential equation 
8(8—3p,—1)(6—3p,—2)u = xu (30) 
formed by extracting those parts of (19) which involve x only, and 
— u = Ai(pi)AZ(p2)w. 
This substitution (in 2) corresponds to the substitution (21) in 2, y, 
and we similarly find that u satisfies (30) if w is any solution of 
8(8—1)(5—2)w = aw, 
dw 


i.e. of i = w. (31) 
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Taking the solution w = e*, we get for the equation (30) the solution 
u = Al(p,)AZ(po)e”, 
which is exactly (x). We similarly find that »(y) is a solution of 
the ordinary equation 
5’ (8’—3q,—1)(8’ —3q.—2)w’ = y*w’, (32) 

extracted from (19) by suppressing 2, 5; and that (x), 4(y) are 
solutions of the equations adjoint, in the operators 5, 5’, to (30), (32), 
namely 3(8+3p,+1)(8-+3p_+2)v = (—a)*v | 

8’(8’+-3q,+1)(8’+-39,+2)v = (—y)v J’ 
which can, of course, be formed out of the original equation (18) 
adjoint to (19). 

Conversely, with &, n, &, 7 so defined as solutions of their respective 
differential equations, we see readily that (29) satisfies the conditions 
(A) postulated for the fundamental solutions and that this is true 
even when 7), Po, 71; J. ave no longer positive integers. More generally, 
if E(x), n(y) are solutions of the ordinary equations 


(33) 


f(—d)E = a, go (—8’)nQ = 9, (34) 


where f, g are polynomials of equal degree m, and &(x), 7(y) solutions 


of their adjoint equations 
SSE =a, 98’) = y™, (35) 
then it is clear that 


" (+) 


| E(at)n(wyt)E(Xt)7(wY te" dt (36) 


satisfies the conditions (A) for fundamental solutions of the partial 


differential equation (5). This is analogous to what has been proved 
in (3) for equations of the form (1). 


7. The null condition for the contour-integral solutions 

So far, then, to satisfy the condition (A), it has been sufficient 
that £, », &, 7 should be any solutions of the differential equations 
(34) and their adjoints. To satisfy the null condition (B) it is neces- 
sary. to be more precise and to link suitably the y, 7 that appear in 
the set of fundamental solutions. 

Returning to third-order equations, let us replace the particular 
solutions 7(y), n(wy), 7(w?y) of (32) by any linearly independent set 
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11> Ne Ng. With these we associate, as is usually the case, the parti- 
cular adjoint solutions* 7, 72, 7; obtained by dividing the cofactors 


(8m, Sy. S73 
a Te 
of the Wronskian 
W =| 8y, 8°, 8y 5 
| 8, 82 43 
m 13 
by the Wronskian itself. Then the sum } 7(wyt)7(w¥t) which 
appears in > U,, will be replaced by the expression 
W) m(yt)  ne(yt) —nalyt) 
6,(¥t) Sy(¥t) Sn3(¥t) , (37) 
m(Yt) (Yt) (Yt) 
which is evidently invariant when 7, 7, 73 are replaced by any 
linear combination of them (of course, with constant coefficients). 
In other words (37) is independent of the particular set of solutions 
"i> Ne» Ng, and we may therefore replace them by any other set with- 
out affecting > U,,. Now the theory of solving hypergeometric 
ordinary differential equations in power-series shows us that the 
equation (32) has, in general,+ three solutions of the form 
Py’), yO Ply®), yy B(y*), 
where each P(y*) denotes a power-series of the form 
ay +a, y®+a, y®+..., 
i.e. a power-series composed of non-negative integer powers of y’. 
If we substitute these three solutions for ,, 72, 73 in (37), we see 
that, regarded as a function of ¢, it is also of the form P(é). Hence 


now > U,, has the form 
(0+) 


. . | &(at)E(Xt) P(e)t? dt. (38) 
6772 


By a second appeal to the theory of solution in power-series we 
see that ¢(at), (Xt), as functions of ¢, will have the forms 


P,(8)+8?:+1P, (8) + 8?2+2P, (8), 
P, (#3) +-t-GP DP, (8) + 1-GP2+2) P(e). 


* Cf. (6), 382 (25). 

+ That is to say, unless 3g,+1, 3q,+2 or their difference 3(¢,—gq,)—1 is 
zero or a positive or negative integer of the form 3n, in which case logarithmic 
solutions appear. 
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Hence &(xt)&( Xt) P(é) will contain no term in t-* unless, exceptionally, 
either 3p,+1, 3p,+2 or their difference 3(p,—p,)—1 is zero or a 
positive or negative integer of the form 3n. Thus, in general, the 
integrand of (38) has zero residue at t = 0, and so > U,, = 0, which 
is the null condition. 

In the exceptional cases in which either (or both) of the differential 
equations (30), (32) has a logarithmic singularity, I shall be content 
to suppose, by the principle of continuity, that the null condition 
persists: I have no doubt that a closer analysis would confirm this. 

It should be observed that the foregoing argument about the null 
condition does not require us to suppose p,, Po, Y1, Y2 positive integers, 
and that it evidently extends to equations of any order. 


8. The asymptotic values of (x), »(y) 

To determine the behaviour of the contour-integral for U,, along 
the associated characteristic it is (naturally) necessary to fix the 
solutions €(x), n(y) of the differential equations precisely, except for 
possible constant multipliers which can be adjusted in the final pro- 
duct £7. I shall define these solutions by means of their asymptotic 
values for large complex x, y. When the £, 7 have been found, we 
can then determine the ‘associated’ £, 7 described above. 

As we have already seen, when p,, p, are positive integers, an 
appropriate solution € is 

(x) = A,(p)A,(pe)e"; 
where the operators A are of degrees p,, p, in 6 and have leading 
coefficients unity. Thus their product can be expanded in the form 
A.(pyA,(ps) ae arte t-C, toot orient. AO, 
where the C’s are constant. Thus 
&(x) --: e7(xPitP2+ OC} gertPet 4... $C, sy) 
= xPitPret P(x), 

where P, denotes a polynomial function of its argument with leading 
coefficient unity. Thus, if x is complex and |x| + 00, 

E(x) rH xP it Pre, (39) 


When 7,, ~, are not both integers, the differential equation (30) is 
no longer soluble in finite terms, the polynomial P,(x~") being re- 
placed by an infinite series: and, what is worse, this infinite series 
is divergent. Fortunately, the asymptotic form (39) remains and is 
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sufficient to identify the precise solution. We could show this by 
extending to equations of higher order a method given by E. L. 
Ince* to determine the asymptotic form of solutions of analogous 
(but more general) equations of the second order. I shall be content 
here to notice that, if in (30) we put 
= Pit Per 7, 

the equation can be written 

(p> “. 


| 
— D—— ei = 0, 
x* x3 


where 6, c are constants, i.e. 


(D?+.3D+3)Da = 1(6p+0+°)a. 
a? x 


The application of Cauchy’s existence-theorem to an equation of 
this form indicates the existence of a solution which converges to 
an arbitrary constant, say unity, as |2| > oo. 

Accordingly, I define, in terms of asymptotic form, the solutions 


. (40) 


é, —_ Pit+Prer é. pr LPitPrgwx é, eae yPitPrgw*x 
Mm ~ yBite, — ng ye, ny ww yttte'y 
Substituting these asymptotic forms in the Wronskian of £,, &, &3 
we at once have for the ‘associated’ adjoint solutions 
Ew fa rtortte-2 ES fear tr te—wr, Ey we fw2rlPrtrtiie-w'sx, 
(41) 
with corresponding asymptotic forms for 7,, 72, 73-T 
Now take the set of fundamental solutions in the modified formt{ 
(0+) 

U, = 90; | Alen XOH (Oe de (r =1,2,8), (42) 
where the characteristic (x —X)+w(y—Y) = 0 now corresponds to 
the suffix 2. To determine the value of U, along this characteristic 
I expand the contour of integration into the ‘infinite’ circle: in virtue 
of the asymptotic forms (40), (41) we then have 


U, ~ = | (art)P1+P2(yt)%+42( Xt) Pr t+2+2)(Y ft) ar tart ele twy-X —w¥ yp? df, 

‘a sare 

* (7), 169-71 (§ 7.31). 

+ These asymptotic forms agree with those derivable from the expressions 
given for é, 7 in § 6. 

t In particular a factor } has been removed since it now appears in the 
asymptotic forms (41). 








234 T. W. CHAUNDY 
Define Wy = e-@rt dy “Gita XPitP2t2Vaitdet2 (7, 


w —— 
so that W~ | pees -arF 


67t e’ 


ow, w elias it 
t 


and therefore - -- —. 
Ox 670 


Hence, along the characteristic x+wy = X+oY, 


and so, by (43), 
6U, 


Cx 


U, — 0, —_ Swarr tPryh +42 X¥ —(1+P2t+2)Y “i +d2+2), 

These (with the similar result for the y-derivative of U,) are the 

‘characteristic’ conditions (C) for the equation we are considering.* 
It is not difficult to see that this argument extends to the corre- 

sponding equation of any order. 


9. An explicit form for the functions £(x), etc. 

Explicit expressions can be obtained for the functions €, », é, 7 in 
the form of (m—1)-fold integrals. I consider only functions of the 
third order and, in particular, equation (30) 

5(5—3p,—1)(6—3p,—2)u = xu, (30) 


where I suppose the factors on the left arranged in descending 


order, i.e. 3p,+2 > 3p,+1> 0. (44) 
In the first instance restrict p,, p. to be integers and operate on the 
differential equation with 

x~8P2+2)(§ —3)...(8—3p,—3) . (6B—3p,—4)...(8—3p,—4). 


This gives, after reduction, 
{8(85—1)(8—2)—a3}a-Br2t?) x 
x 6...(8B—3p,). (6—3p,—1)...(6—3p.—1l)u = 0, 
where the operators in each group descend by differences of 3. Thus, 
for a suitably chosen solution u = £, we have 
x-BP2+2§_, (S—3p.).(6B—3p,—1)...(8—3p,—1)é = e*, 
d — d Pe Pt" ant = idl 
x*dx x*dx = 


* Cf. equation (28) above. 
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Now it is easily verified that a particular solution of the differential 


equation a: 
(sre) fe) = 9) 


x*dx 


3 n! ’ 


* 3___ 43\ n #2 
is given by f(z) = [ af? *) t? dt 
h 


where h is any convenient constant. 
If we take h = —oo and apply this formula twice in the differential 
equation (45), we get at length 


x t 
3 _— f3\P2—-P1  $-(3P2+2) 3__3\ De 
ear) = gdp [ (" “) t | [ e( “| rdrdt 
F. 3 (P2—Pr)! J 3 Pp! 


We suppose the real part of 2 negative to avoid difficulties at ¢ = 0 
on the path of integration in ¢. For more general values of 2 this 
double integral would need to be extended into a double contour- 
integral of the type introduced by Hankel for Bessel’s functions. 

Removing the integer-restriction on p,, p, but retaining the in- 
equalities (44), we can replace this double integral by 


o(: x t 
~3p,+1 e 3__ #3\ De—P1 3__ -3\ De 
= — = (" *) ¢-Gp2+2) [ e( =| tdrdt. 
D(ps—py+1)P(p.+1) , 3 3 
— (47) 


So far I have proved (47) a solution of the equation (45) but not 
necessarily of the simpler original equation (30). To complete the 
argument write 


(46) 


. 


t= 2/t'r’ 
and drop accents. Then, after reduction and extruding a certain 


constant C, we have 
1 2 
& — (y3p2+2 [_[ enteaa—t-syprnie-t(1—2- sp 


. 


0 0 
Y 2t+2¢" 29 17 
= Crsrst2t’ gay, 


dtd 


tr 


11 
where v= | [ emoow 


. 


0 0 
and f(t) = (Lt), f(r) = (1-79), 
We can easily verify that ¢, % satisfy the respective equations 
(5+ 3p2—3p,+1)h = (8,—2)84, (8, + 3p. +2) = (8,—1)7*Y. 
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Then é’ is a double integral of Mellin’s type such that, if certain 
expressions vanish at the limits 1, 0, 
8(8-+3p,—3p,+1)(6+3p,+ 2)é’ 


— 3(8—1)(8—2) [ [ ee (Oyr8f(x) Oe = xf’. 


~ « 


0 0 

Thus, from (48), 6(8—3p,—1)(—3p,—2)é = aE, 
if the limits in the double integral have been well chosen. 

It is perhaps easier to interpret the integral for £’ as a double 
Laplace integral, writing the differential equation for €’ as 

Ag’ = {a?(D§—1)+3(p+q+2)aD?-+ (3p+3)(3q¢+2)D}é’ = 0, (49) 
where p= Po, ¢ = Po—P;. Writing also ¢t = 1/t’, 
integral for €’ and dropping accents we have 


>= 1/7’ in the 


é* = [ [ e%(1—#)(1—15)Pr dtdr. 
Substituting this in (49) we obtain, on reduction, 


8 


Ad’ = —x | (1—73)P[ew(1—88)2+1]° dr — 
1 


~ | t(1—18)¢[e""(1—73)P"(atr-+3q+2)]” | di, 
1 


which vanishes when x < 0, p+1, ¢+1 > 0, the last two inequalities 
being secured by (44). 
To investigate the asymptotic form write in (47) 
7T=t—T, t= xt’ 
and drop accents. Then, after reduction, 


0 « 


é UPit+Prer | ‘ ( ne {8 \P2-P1 
= . = e~ eis —. eer x 
D(po—py+1)0(p.+1) . x =) 


, 7 73 \"" _t —(H2+2) _ ++# 
as if 2—t* 3(@—0? . x : x ac 


Thus, for x of large absolute value, 

© re) 
UPitPer 
é ~ e-l{P2-P1 dt e7T7P: dr = “PitPer, 
0 


D(pe—Py +1) (p_+1) . 
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By (41) this form identifies € as the required £,. We can evidently 
write down analogous double integrals having the asymptotic form 
of €,, €, and therefore to be identified with them.* Further, we can 
write down similar double integrals which are solutions of the 
differential equation for — and have the same asymptotic forms as 
E,, &,, &3; these, then, must be the same as the ‘associated’ ,, &, &; 
as calculated from the Wronskian of £,, &, €;. We deal similarly 
with 7, 7. 


10. The second-order equation 

I conclude with a brief account of the equation of the second order 
for which the theory is somewhat simpler. In the first place, the set 
of fundamental solutions reduces to the single fundamental solution 
of Riemann’s theory: more precisely, there is a pair of such solutions, 
which, in virtue of the null condition, are equal and opposite, and, 
naturally, it is sufficient to pick out one of these and regard it as 
the only fundamental solution U(X, Y; x, y). 

I shall take the partial differential equation and its adjoint to be 


y°d(8+-2p+1)V = 28'(8’+29+1)V, (50) 
y?3(8—2p—1)U = 2°8'(8’—2q—1)U. (51) 
The characteristics are the lines xy = constant, and the charac- 


teristic condition (C) now requires that, along each of the two 
characteristics x—X = +(y—Y), 


U= ePytX —P+Dy a+, (52) 


As before, when p, g are positive integers, we can give the funda- 
mental solution in the form 


f ‘ 7)\2 
/ laa ; (@@—X)4(y—Y)r+% 
U = X-ewey eA, (p)A,(g)Ax(P)Ar(@) oat 


> 


(53) 


— A,(p) = (8—1)(8—8)...(8—2p+1), 
A,(q) = (8’—1)(8’—3)...(8°—-2g+ 1), 


and Ax(p), Ay(q) are similar operators in X, Y. In the operand 


* For no C, €,+C, &,+C; & has, for all 2, the same asymptotic form as an 
&, except, of course, &, itself. 
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of (53) the choice of either sign in + leads to the same expression 
for U: this may be regarded as an aspect of the null condition. 

In expressing U as a contour-integral the differential equations 
defining £, » are 


5(5—2p—1)é = 2°, 8'(8’—2qg—1)n = y?n, (54) 


which are evidently variants of Bessel’s equation. The first of these 
equations has the solutions 


xPttK .4(+2). 


More exactly, we remark that, when p is a positive integer,* 
(5—1)(8—3)...(6—2p+l)e-* = (—)# I ars *K,, 43(2). 

Replacing p+4, g+4, for convenience, by p, v, we can define 
U = =(x), (} ) {x (—at)K,(—yt)K,,(Xt)K, (Yee dt. 


Changing signs in the argument of K,, consider alternatively 


; (0+) 
o = = (5) (§) | Ku—aoKy t)K,(Xt)K,(—Yet)tdt. (56) 


Then U’—U 


+) 

2% y emente . . . > iar 

=3 (5) (¥) { K (—at)K (Xt){K,(yt)K,(—Yt)—K,(—yl)K,(Yo)}¢ dt. 
(57) 

Now on substituting first} for K,(—Yt), K,(—yt) and then{ for 

K, (yt), K,(Yt), we reduce 
K,(yt)K,(—Yt)—K,(—yt)K,(¥0) 

to sin® cosec va{ I, (yt)I_,(Yt)—L(Yt)L,(yt)}, 
and this can be expanded in a power-series P(t*) containing only 
non-negative integer powers of ¢*, since each term can be separately 
so expanded.§ K,(—at), K,(Xt) can each be expanded in power- 
series of the form é#P,(é?)+?¢-“P,(#), and so, as in the argument of 
$7, the integrand of (57), in general, has zero residue at ¢ = 0, and 
accordingly U = U’. This equivalence, in the second-order equation, 


* Cf. (8), 79 (6) and 80 (13). + From (8), 80 a 
t From (8), 78 (6). § Cf. (8), 77 (2). 
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of the two forms (55), (56) for the fundamental solution, corresponds, 
of course, to the null condition in equations of higher order. 
Expanding the contour in (55) into the ‘infinite’ circle and so 
replacing the Bessel’s functions by their asymptotic forms,* we 
obtain (returning from p, v to p+4, q+4) 
dt 


Uw~ . myiX-osny ern | deren’, (58) 
274 t 


Thus, along the characteristic <+y = X+Y, 
a sue aPytX —v+by —a+); (59) 


similarly U’ has the same value along the other characteristic 
x—y = X—Y. In view of the equivalence of U, U’ this establishes 
the ‘characteristic’ condition (C).t 

For purposes of reference it is perhaps better to use a substitution 
of the type (174), (178) to transform the equation (51) into a 


‘normal form’ y?(8°—p2)U, = 2(8"2—v2)U, (60) 


with the fundamental solution 


(0+) 
3 
U=5 | K,(—2t)K,(Fy)K,(X)K,(+Yot dt. (61) 


By a substitution of the same type we can transform the equation 
mio 6+ p)(8—P—1]U, = 2*(8'+-q)(8'—q—1) Uy 
_ 80, _ fete _ eat Diy, wb 


= y 
The change of independent variables 
L=H4tY,p YrurN—-n 


n2qT 
then gives 4 Uy = (eee 1) -§ | U;. (63) 
02, OY, (+H)? (%— y:)") 

* From (8), 202, § 7.23 (1). 

+ Reference to the condition as formulated in § 2 will suggest that a factor 
4 has been omitted in (59). Actually this is accounted for by a slight difference 
between the technique that I use for the equation with m characteristics and 
Riemann’s method applicable to the equation with two characteristics. To 
reconcile them we could here have taken the fundamental solution to be 


3(U+U’). 
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Now this is a case of Darboux’s equation which I have discussed in 
(5), giving its fundamental solution (in the notation of that article) as 
m= Al? 4 
2 2 é ” 

where £= — (%—X)(Hi—N) n= (%—X,)ii—Vi) 
(%,+%)(X1 +H) (%,—y1)(X,—N) 
Since the fundamental solution of a partial differential equation of 
the form s = Az is known to be unique, we can identify this expres- 
sion for U, with that obtained from (55), thus getting (in the more 
convenient variables x, y)* 
(0+) 
2% ry) [ x 4 r aK (Y 1 |P @ 
SNOUXY) | Kya — 20K ya(—Y Ky KOK (WOM = ALP 2 
(64) 

Xf—ig—FF aaa (z—X)*—(y—Y)? 

4x.X ; 4yY : 


By using one of the expressions for /\, given in (5), for example 


where €= ite (65) 


p 
Al 
"le 9 
we express (64) as an identity connecting Bessel’s and Legendre’s 
function (which I have not so far found in the literature). 


J si (es +9541) | P,(1—2€-+ 2¢t)P,(1—2nft) dt 
0 


* In transforming from U to U, we have also to make an adjustment for 
the change from 4, 5’ to @/éx, é/éy as the operators defining ‘adjoint’. The 
final form can be checked by noticing that the equations (62), (63) are self- 
adjoint and that accordingly U, should be symmetrical in (x, y) and (X, Y). 

t This is (5), 237 (9). 
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